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ABSTRACT 


A  two  layer  ocean  mixed  layer  model  is  used  to  study  the  relationship  between  the 
onset  of  freezing  and  mixed  layer  depth,  wind  forcing,  surface  buoyancy  flux,  and  tem¬ 
perature  and  salinity  changes  between  the  two  layers.  Universal  non-dimensional  pa¬ 
rameters  for  stability  and  surface  forcing  are  derived  and  related  to  the  maximum 
freezing  rate.  A™ lytic  solution:  to  the  model  are  found  in  terms  of  tne  universal  pa¬ 
rameters  and  the  model  turbulent  mixing  tuning  constants. 

Sensitivity  studies  show  the  dependence  of  the  freezing  rate  on  the  stability  as  de¬ 
fined  by  the  salinity  and  temperature  jumps,  the  forcing  by  the  wind  stress  and  the  sur¬ 
face  heat  flux,  and  the  mixed  layer  depth.  Results  show  that  an  increase  in  heat  flux 
produces  a  nearly  linear  increase  in  the  freezing  rate.  Mixing  energy  from  the  wind, 
proportional  to  the  wind  speed  cubed,  results  in  a  nearly  linear  decrease  in  freezing  rate. 
There  is  a  non-linear  relationship  between  the  temperature  jump  between  layers  and  the 
freezing  rate.  A  warming  of  the  deeper  layer  decreases  the  freezing  rate  and  ultimately 
prevents  freezing.  A  non-linear  relation  was  also  found  between  the  salinity  jump  and 
the  freezing  rate.  An  increase  in  the  deeper  layer  salinity  causes  an  increase  in  the 
freezing  rate  and  leads  to  the  maximum  expected  freezing  rate.  A  nearly  hyperbolic  re¬ 
lationship  between  the  mixed  layer  depth  and  the  freezing  rate  was  found.  As  the  mixed 
layer  depth  increases  from  near  zero,  there  is  a  rapid  increase  in  the  freezing  rate,  and 
then  the  maximum  expected  freezing  rate  is  approached  asymptotically. 

A  relationship  between  the  forcing  parameter  and  the  stability  parameter  was  dp- 
rived  which  defined  areas  where  freezing  could  occur.  For  low  values  of  the  stability 
parameter  there  is  the  possibility  of  freezing  for  any  value  of  the  forcing  parameter.  For 
larger  values  of  stability  there  can  only  be  freezing  for  small  values  of  the  forcing  pa¬ 
rameter.  There  is  a  critical  value  of  the  stability  parameter  above  which  freezing  is  not 
possible  regardless  of  surface  forcing. 
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I.  INTRODUCTION 


The  production  of  sea  ice  is  significant  on  most  scales  of  interactions  between  the 
Arctic  oceanic  and  atmospheric  boundary  layers.  The  Arctic  ice  pack  and  marginal  ice 
zone  (MIZ)  have  been  studied  to  determine  on  various  scales  the  interactions  of  the  ice 
on  the  local  environment  and  the  importance  to  the  global  climate.  These  studies  range 
from  ocean  wave-ice  interactions  (Martin  and  Kauffman,  1981,  Christensen,  1983, 
Squire,  1984,  Wadhams  et  al.,  19SS).  internal  wave  activity  (Muench  et  al..  1983, 
Morison,  1985,  Levine  and  Paulson.  1987,  Sandven  and  Johannessen.  1987,  Foster  and 
Eckert,  19S7,  McPhee  and  Kantha,  19S9,  Eckert  and  Foster,  1990),  inertial  oscillations 
(McPhee,  1978,  Lagerloef  and  Muench,  1987),  ice  edge  banding  (Martin  et  al.,  1983, 
Wadhams,  1983,  Smedstad  and  Roed,  1985,  Chu,  19S7),  leads  (Andreas,  19S1,  Kozo, 
1983,  Chu,  19S6,  Gow  et  al.,  1990),  polynyas  (Topham  et  al..  19S3,  Smith  et  al.,  19S3. 
Flartog  et  al.,  1983,  Chu,  1986),  ice  floe  collisions  (Shen  et  al.,  19S7,  Martin  and  Becker, 
19S7,  19SS,  Lu  et  al.,  1989),  sea  ice  motion  (Thorndike  and  Colony,  19S2,  Colony  and 
Thorndike,  1984,  1985,  Thorndike,  1987,  Thompson  et  al.,  19SS,  Scrreze  et  al.,  1 9 S 9 ) , 
and  basin  wide  climate  studies  (Washington  et  al.,  1980,  Hibler  and  Walsh,  1 9S2. 
Niebauer,  1983,  Semtner,  19S7,  Nakamura  and  Oort,  1988,  Ingram  et  al.,  1989, 
tiakkinen  and  Mellor,  1990,  Aagaard  and  Carmack,  1990). 

On  smaller  scales,  it  has  been  found  that  for  newly  formed  sea  ice,  frazil  and  grease 
ice,  the  drag  coefficient  can  be  decreased  to  one  half  of  that  of  the  open  ocean  (Macklin, 
1983,  Guest  and  Davidson,  1991).  On  the  other  hand,  a  build  up  in  thicknec<:  and  in¬ 
crease  in  surface  roughness  from  lateral  stresses  can  increase  the  drag  coefficient  from  2 
to  3  times  that  of  the  open  ocean  (Macklin,  1985,  Pease  et  al.,  1983,  Walter  et  al.,  1984, 
Overland,  1985,  Guest  and  Davidson,  1987,  Guest  and  Davidson,  1991).  The  effect  on 
wind-stress  imparted  momentum  to  the  ocean  and  ice  can  dominate  the  flow  of  the  sea 
ice  and  the  underlying  surface  water. 

In  higher  latitudes  the  ocean  moderates  the  atmospheric  temperature,  with  air  tem¬ 
perature  being  close  to  sea  surface  temperature.  While  the  ocean  is  ice  free  there  can 
be  a  large  heat  flux  into  the  atmosphere.  As  sea  ice  forms  and  decreases  the  amount  of 
open  water  exposed  to  the  atmosphere,  the  heat  flux  to  the  atmosphere  will  decrease 
significantly.  There  is  a  similar  alfect  on  the  moisture  flux  from  the  sea  surface  to  the 
atmosphere  as  the  areal  coverage  of  sea  ice  extent  increases.  This  results  in  seasonal 
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changes  in  clouds  and  fog  over  the  ice,  with  a  maximum  in  cloud  cover  in  summer  and 
a  minimum  in  winter  (Sater  et  al..  1971). 

The  formation  of  sea  ice  also  influences  the  radiation  budget  of  the  atmosphere. 
The  differences  in  albedo  between  open  water,  thin  ice,  thick  ice,  and  snow  covered  ice, 
will  change  significantly  the  penetration  of  incoming  solar  radiation  into  the  ocean  as 
well  as  reradiation  back  to  the  atmosphere  (Grenfell  and  Perovich.  19S4,  Maykut,  19S6, 
Ross  and  Walsh.  1987,  Ingram  et  al.,  19S9).  Longwave  energy,  coupled  with  the 
moisture  of  the  atmospheric  boundary  layer,  is  greatly  influenced  when  open  water  be- 
comes  covered  with  sea  ice  and  snow. 

In  the  past  decade  there  has  been  increased  investigation  into  the  formation,  extent, 
and  compaction  of  sea  ice.  Ice  physics  studies  have  shed  new  light  on  the  mechanisms 
of  ice  formation  (Hanley,  197S,  Omstedt  and  Svensson,  19S4),  growth  (Bauer  and 
Martin,  1983.  Wakatsuchi  and  Ono,  1983).  and  compaction  (Lepparanta  and  Hibler. 
1985)  on  many  scales.  Modelers  have  progressed  from  one  dimensional  models  to  three 
dimensional  dynamic  and  thermodynamic  models  of  varying  scales  of  the  ice  pack  and 
the  MIZ  (Pollard  et  al.,  19S3,  Mcphee,  1983,  Roed,  19S4.  lkeda,  19S6,  Mellor  et  al., 
1986,  McPhcc.  19S7,  MePhee  et  al..  19S7,  Semtner.  19S7.  Houssais.  1988,  Kamha  and 
Mellor,  1989,  Mellor  and  Kantha,  19S9,  lkeda,  19S9,  Lu  et  al.,  1990).  These  model 
studies  have  improved  the  understanding  of  the  formation  of  deep  water  from  freezing 
(Carmack  and  Aagaard.  1973.  Swift  and  Aagaard,  19S1,  Hakkinen.  19S7.  Swift  and 
Koltermann,  1988),  salt  chimneys  (Killworth.  1979).  eddy  formation  and  interaction 
(Manley  and  Hunkins.  19S5,  Hakkinen.  19S6,  Johannessen  et  al.  19S7,  Smith  and  Bird, 
19S8,  Padman  et  al..  1990),  and  other  small  to  large  scale  features  found  in  the  arctic. 
Unfortunately,  few  studies  have  been  done  on  the  effects  of  the  dynamical  properties  of 
the  mixed  layer  and  their  various  interactions  which  contribute  to  ice  formation,  or 
prevent  ice  formation.  Chu  and  Garwood  ( 19S8)  have  indicated  that  there  may  be  some 
coupled  ice-ocean  domains  where  the  effect  of  entrainment  and  surface  buoyancy  flux 
may  not  be  fully  understood. 

This  study  will  focus  on  various  thermodynamic  and  dynamic  properties  of  the 
mixed  layer,  wind  stress  and  buoyancy  effects  to  better  understand  the  air-sea  coupling 
leading  to  the  formation  of  sea  ice.  The  identification  of  regimes  where  freezing  can 
occur  will  be  of  special  interest. 

Chapter  II  introduces  the  theory  of  mixed  layer  dynamics  and  thermodynamics  in¬ 
volved  in  the  freezing  of  salt  water.  Chapter  III  presents  a  simple  one-dimensional 
model  with  associated  thermodynamic  equations.  Chapter  IV  shows  the  interaction  of 


the  various  dynamic  and  thermodynamic  processes,  and  solutions  are  developed  for  the 
problem  of  finding  regimes  of  freezing. 
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II.  THEORY 


A.  MIXED  LAYER  DYNAMICS 
1.  The  Oceanic  Mixed  Layer 

The  oceanic  mixed  layer  is  a  quasi-homogeneous  layer  in  the  ocean,  from  the 
ocean  surface  to  a  depth  called  the  mixed  layer  depth,  in  which  the  various  properties 
of  the  -water,  thermal  energy,  and  salinity  are  well  mixed  by  the  turbulence. 


Temperature  ('Q  5slinll>  {7£  ) 

S.O  7 . 5  10. 0  12. S  1S.0  17.5  20.0  22.5  25.0  27.5  30.0  32.5  35.0 


Figure  I.  Sample  Midlatitude  Temperature,  Salinity,  Depth  Profiles:  The 
temperature  profile  shows  a  uniform  Mixed  Layer  temperature  and 
decreasing  temperature  with  depth.  The  salinity  profile  shows  a  uniform 
mixed  layer  salinity  and  a  gradual  increase  with  depth.  (Jessen  et  al., 
1989) 


A 


Pressure  (dbar) 


Figure  2.  Typical  Arctic  Temperature  And  Salinity  Depth  Profile:  The 
temperature  profile  shows  a  cold  shallow  mixed  layer  over  a  warmer 
deeper  layer.  The  salinity  profile  shows  a  less  saline  mixed  layer  over  a 
more  saline  deeper  layer. 
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Figure  17.  VERTICAL  DISTRIBUTIONS  OF  TEMPERATURE 
AND  SALINITY  AT  SIX  LOCALITIES 


Figure  3.  Vertical  Distribution  of  Temperature  and  Salinity  at  Six 
Locations:  Salinity  and  Temperature  profiles  at  six  locations  in  the 
Arctic  Ocean  and  adjacent  seas.  All  temperature  profiles  show  a  cooler 
mixed  layer  overlying  a  wanner  layer,  and  all  salinity  profiles  show  less 
saline  mixed  layer  water  overlying  more  saline  water.  (Courtesy  Sater  et 
al.,  1971) 
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Figure  4.  MIZEX  87  Temperature,  Salinity,  Density  Profile:  Temperature, 
Salinity  and  Density  profiles  from  R/V  HAKON  MOSBY  taken  at  78.1 
°N  2.5°W  on  3  April  1987.  The  temperature  profile  shows  a  uniform 
mixed  layer  temperature  of  -1.8  °C  to  80  meters  and  a  rapid  increase  to 
1.0°  C  at  120  meters  and  then  decrease  with  depth.  The  Salinity  profile 
shows  a  relatively  uniform  salinity  of  34.5  in  the  mixed  layer  increasing 
to  34.95  at  150  meters  and  maintaining  a  relatively  constant  value  below 
150  meters.  The  Density  profile  shows  a  large  density  increase  just 
below  the  mixed  layer.  (Sandven,  et  al.,  1987) 
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Over  much  of  the  ocean  the  mixed  layer  temperature  profile  is  similar  to  Figure 
1,  taken  off  the  California  coast  (Jesscn  et  al.,  1989),  where  the  mixed  layer  is  warmer 
than  the  underlying  water.  In  the  mid-latitudes  where  the  amount  of  evaporation  is 
exceeded  by  precipitation,  the  mixed  layer  salinity  is  lower  than  that  of  the  underlying 
water  column,  as  is  evident  ir  .  igure  1.  Here  the  water  density  difference  between  the 
mixed  layer  and  the  lower  layer  is  a  function  of  both  temperature  and  salinity.  However, 
in  the  Arctic,  the  mixed  layer  is  typically  colder  than  the  underlying  water  as  shown  in 
Figure  2.  For  the  arctic,  the  salinity  of  the  mixed  layer  is  typically  less  than  that  of  *he 
underlying  water  as  shown  in  Figure  2.  (Figure  2  show's  hypothetical  conditions  that 
are  simplified  for  use  in  a  one-dimensional  mixed  layer  model.) 

For  temperatures  close  to  the  freezing  point,  changes  in  water  density  are  mostly 
a  function  of  salinity.  This  allows  a  cooler  less  saline  mixed  layer  to  be  much  more 
buoyant  than  a  warmer  more  saline  layer.  Therefore  salinity  dominates  the  density 
gradient,  allowing  colder  water  to  overlie  warmer  water  (cooler  water  is  normally  denser 
than  warmer  water).  Depending  on  the  season  and  the  proximity  to  fresh  water  from 
river  run  ofT,  the  mixed  layer  may  be  near  freezing  and  up  to  4  °C  cooler  than  the 

g 

underlvine  water.  The  salinitv  also  mav  be  up  to  5  - —  less  in  the  mixed  laver  than  in 

'  '  kg 

the  layer  below.  Figure  3  shows  a  variety  of  temperature  and  salinity  profiles  from  six 
different  locations  in  the  Arctic  Ocean  and  outlying  seas.  The  general  conditions  of 
cooler  less  saline  waters  overlying  warmer  more  saline  waters  appears  in  all  six  of  the 
regions  depicted  in  Figure  3. 

The  energy  needed  to  mix  the  surface  layer  properties  comes  from  two  sources. 
The  normally  dominant  source  is  the  stress  imparted  by  the  wind  to  the  water  or  ice 
surface.  As  the  wind  speed  increases  the  energy  imparted  to  the  mixing  process 
increases.  The  second  source  of  energy  is  due  to  density,  or  buoyant  instabilities  in  the 
water  column.  Buoyant  instabilities  are  caused  by  three  sources:  heating  from  below, 
cooling  from  above,  and  downward  salinity  flux  at  the  surface  or  at  the  ice-water 
interface.  The  salinity  fluxes  occur  when  there  are  net  changes  in  the  mixed  layer  due 
to  freezing,  melting,  evaporation,  or  precipitation  at  the  surface. 


8 


2.  Freezing  Processes 


In  the  Arctic  and  Antarctic  waters,  a  typical  temperature,  salinity,  and  density 
profile  would  show  cooler,  less  saline  and  less  dense  water  overlying  wanner,  more  saline 
and  more  dense  water.  Figure  4  shows  temperature,  salinity  and  density  profiles 
obtained  from  the  R/V  HAKON  MOSBY  during  the  MIZEX  87  experiment  (Sandven 
et  al.,  1987),  demonstrating  these  properties.  These  profiles  will  be  used  later  as  a 
guideline  for  model  initial  conditions  in  chapter  4.  During  the  winter  freezing  period, 
very  cold  air,  — 10  °  C  to  —25  °C,  usually  overlies  the  ice  and  the  open  ocean  near  the  ice 
edge.  This  causes  a  heat  flux  from  the  ocean  surface  to  the  air.  In  the  case  where  ice 
is  already  present,  the  heat  flux  is  from  the  ocean  through  the  ice  to  the  air.  This 
upward  heat  flux  cools  the  mixed  layer  to  the  freezing  point  [where  the  freezing  point 
is  mainly  a  function  of  salinity].  Also,  heat  may  flux  from  the  warmer  underlying  waters 
into  the  mixed  layer,  which  may  slow  down  the  cooling  of  the  mixed  layer,  or,  depending 
on  the  amount  of  heat  flux  from  below,  it  may  even  warm  the  mixed  layer. 

Bauer  and  Martin  (1983)  describe  two  conditions  when  freezing  occurs  in  open 
leads  or  open  ocean,  for  both  windy  or  calm  conditions.  Under  calm  conditions  a  thin 
horizontal  sheet  of  ice  forms  and  grows  downward.  Under  these  conditions  the  major 
thermodynamic  process  is  long-wave  radiative  cooling  of  the  surface  leading  to  cooling 
and  freezing. 

Under  windy  conditions  small  disc-shaped  crystals,  called  frazil  ice,  are  formed 
and  driven  down  wind  where  they  can  become  concer»*rated  and  freeze  together  in  the 
form  of  grease  ice.  In  leads,  Bauer  and  Martin  (1983)  have  shown  that  this  process  of 
ice  formation  will  grow  initially  downward  to  a  minimal  depth  and  then  grow 
horizontally  to  the  upwind  edge  and  eventually  close  off  the  lead.  Under  windy 
conditions  the  major  thermodynamic  processes  for  freezing  are  from  the  latent  heat 
fluxes  and  sensible  heat  fluxes. 

Uewis  and  Perkins  (1983)  describe  conditions  where  freezing  occurs  under 
already  existing  ice.  This  happens  when  the  ice  surface  near  the  ice-ocean  interface  is 
colder  than  the  freezing  point,  and  ice  crystals  form  and  grow  downward.  They  describe 
a  process  where  in  some  cases  the  mixed  layer  undergoes  an  upward  heat  flux  until  a 
condition  called  supercooling  occurs.  For  this  case  the  water  at  the  interface  is  actually 
colder  than  freezing  by  about  .01  °C,  and  only  then  will  freezing  begin.  (Lewis  and 
Perkins,  1983;  Omstedt  and  Svenson,  1984). 
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3.  Heat  Fluxes  in  the  Mixed  Layer 


The  heat  fluxed  across  the  water  surface  or  ice  surface  generates  and  sustains 
the  freezing  process,  and  so  surface  heat  flux  is  the  primary  parameter  considered  in  the 
formation  of  ice.  The  various  mechanisms  that  produce  this  heat  flux  have  different 
effects  on  the  mixed  layer.  The  largest  component  of  the  heat  flux  process  at  an 
ice-ocean  interface  is  the  sensible  heat  flux.  The  primary  heat  lux  processes  at  the 
air-ocean  or  air-ice  interfaces  include: 

•  Shortwave  Radiation 

•  Longwave  Radiation 

•  Latent  Heat  Flux 

•  Sensible  Heat  Flux 

•  Conductive  Heat  Flux 

The  latent  and  sensible  heat  fluxes  are  dependent  on  air-sea  differences  and  on 
wind  stresses  at  the  air-ocean  or  air-ice  interfaces.  The  radiative  heat  fluxes  depend 
more  on  cloud  cover,  i.e.  clear  skies,  fog,  low  clouds,  high  clouds,  etc. 

a.  Shortwave  radiation 

Shortwave  radiation,  or  solar  radiation,  is  a  downward  radiation  into  the  ice 
or  the  ocean  mixed  layer.  Mavkut  (1986)  has  parameterized  the  net  shortwave  energy 
flux,  Fr ,  as 

“  *c3Xl  ~  «)  (2.1) 

where 

•  F„  -  incoming  solar  radiation  (function  of  sun  angle) 

•  k  =  seasonal  cloudiness  parameter  ranging  from  .15  in  March  to  .5  in  August 

•  c  =  cloudiness  in  tenths 

•  a  =  albedo  of  the  surface  (function  of  wavelengths) 

Because  both  the  ice  and  ocean  are  translucent,  the  shortwave  energy  can 
penetrate  the  surface.  Mavkut  (1986)  has  shown  a  relation  for  the  heat  flux  Ft  at  depth 
z,  of: 
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(2.2) 


F  =  i  F e~K(2~A) 
where 

•  z  =  depth  (positive  downward) 

•  /„  =  percentage  absorbed  in  a  thin  surface  layer  (function  of  surface  type  (  i.e.  ice 
or  water.  For  ice  and  z  >  .1  m,  i,  is  about  .3),  and  frequency  (more  is  absorbea  in 
the  red  and  infrared  than  blue  and  ultraviolet)) 

•  Ft  —  net  shortwave  heat  flux  at  the  surface 

•  k  =  bulk  extinction  coefficient  of  the  medium 

Clouds,  sun  angle  and  surface  type  will  determine  how  much  energy  is 
incident  at  the  water  or  ice  surface,  and  turbidity  and  ice  type  determine  how  much 
energy  passes  through  to  the  mixed  layer  and  below. 

b.  Longwave  Radiation 

Longwave  radiation  is  the  net  radiation  between  that  given  off  by  the  water 
or  ice  surface  to  the  atmosphere,  and  that  absorbed  by  the  ice  or  water  surface  when 
radiated  by  the  atmosphere.  Maykut  and  Church  (1973)  parameterized  the  downward 
longwave  radiation,  F,d  ,  by: 

F{d  =  .7855(1  +  .2232c2'75)ct7^  (2.3) 

where 

•  c  =  cloudiness  in  tenths 

•  a  -  Stefan-Boltzman  constant 

•  T„  =  air  temperature  ( °k ) 

and  upward  longwave  radiation,  Ft  ,  by: 

Ft.-rf  (2-4) 


where: 

•  t,  -  emissivity  of  the  surface  (.97  for  ice,  water,  melt  ponds,  leads,  and  .99  for 
snow) 

•  T,  —  water/ice  surface  temperature  (°k) 
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The  net  longwave  radiation,  F„  can  be  derived  from  equations  (2.3)  and  (2.4)  as: 
Ft=  .7855(1  +  .2252c1'1 5)oT^  —  .91oT^ 


(2.5) 


With  this  formula,  Maykut  and  Church  found  no  relation  between  the  water 
vapor  and  the  longwave  radiation,  possibly  due  to  the  low  values  of  vapor  pressure 
found  in  the  Arctic.  This  shows  the  role  the  clouds  play,  and  also  the  major  role  of  the 
air  temperature  on  longwave  radiation. 

c.  Latent  Heat  Flux 

Latent  heat  flux  is  caused  by  the  evaporation  of  water  at  the  water  surface 
or  ice  surface.  Maykut  (1986)  has  parameterized  the  latent  heat  flux,  Ft  ,  as 

F  =  pLC’iou-q0)  (2.6) 

where 

•  p  =  air  density 

•  L  =  latent  heat  of  vaporization,  or  sublimation 

•  C,  =  bulk  transfer  coefficient  for  latent  heat  (a  function  of  surface  roughness) 

•  u  =  wind  speed 

•  qa  -  specific  humidity  of  the  air  at  a  reference  level  (usually  10  meters) 

•  q0  =  specific  humidity  of  the  air  at  the  surface 

This  shows  the  major  factors  affecting  the  latent  heat  flux,  the  humidity 
difference  between  the  air  and  the  water  surface,  and  the  wind  speed.  Thus,  latent  heat 
flux  will  have  its  greatest  effect  during  dry  windy  events,  and  its  least  effect  during  moist 
windy  events  or  calm  events. 

d.  Sensible  Heat  Flux 

Sensible  heat  flux  is  the  result  of  turbulent  heat  conduction  between  the 
atmosphere  and  the  water  or  ice  surface.  Maykut  (1986)  has  expressed  the  sensible  heat 
flux,  F, ,  as 

Fs  =  PcpCsu(Ta  —  T0)  (2.7) 
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where 


•  p  =  density  of  the  air 

•  cp  =  specific  heat  of  the  air 

•  C,  —  bulk  transfer  coefficient  for  sensible  heat  (a  function  of  surface  roughness) 

•  u  =  wind  speed 

•  T„  =  air  temperature  at  a  reference  level  (usually  10  meters) 

•  Tc  =  temperature  at  the  surface 

This  shows  the  major  factors  affecting  the  sensible  heat,  the  temperature  difference 
between  the  air  and  the  ice, 'water  surface,  and  the  wind  speed. 

e.  Conductive  Heat  Flux 

For  freezing  to  continue  after  initial  ice  formation  there  must  be  a  heat  flux 
through  the  ice.  This  is  a  conductive  heat  flux,  Fc  ,  and  has  been  parameterized  by 
Maykut  (1986)  as 

F'-H.  4f)„  12-8) 

o  z 


where 

'T* 

•  (  )„  =  the  temperature  gradient  at  the  surface 

•  k,  =  thermal  conductivity  of  the  ice  (a  function  of  temperature  and  salinity  of  the 
ice) 

For  sea  ice,  k,  may  be  expressed  as 
BS, 

kt-{k0  +  A=!-)  (2.9' 

1  i 


where 

•  k,  -  thermal  conductivity  in  pure  ice  (a  function  of  temperature) 

•  p  =  .13  W/m 

•  5,  =  salinity  of  the  ice  in  ppt 

•  T,  -  temperature  of  the  ice  in  °C 

The  effect  of  snow  on  the  ice  is  to  insulate  the  ice  from  the  air  due  to  the 
low  conductivity  of  snow,  which  is  about  one  order  of  magnitude  lower  than  in  pure  ice 
(Maykut,  1986). 
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Figure  5. 


Values  of  a  g  and  f  g  at  Various  Temperatures  and  Salinities:  The 
relationship  between  the  coefficient  of  thermal  expansion,  «*x  10-  ,  at 
different  temperatures  and  salinities  (solid  lines),  and  the  coe  icient 
haline  contraction,  fc  x  10*  ,  at  different  temperatures  and  salinities 

(dashed  lines). 
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Figure  6.  Values  of  a  g  And  P  g  Near  Freezing  Temperatures:  The  relationship 
between  the  coefficient  of  thermal  expansion,  ag  x  104  (solid  lines),  at 
different  temperatures  and  salinities,  and  the  coefficient  of  haline 


contraction,  Pg  x  104  (dashed  lines),  at  different  temperatures  and 
salinities,  near  freezing  temperatures.  At  these  temperatures,  o.g  tends 
to  0,  while  Pg  is  slightly  larger  than  it  is  for  warmer  temperatures. 
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4.  Salinity  Fluxes  in  the  Mixed  Layer 

Salinity  fluxes  in  the  mixed  layer  are  dependent  on  processes  at  the  two 
interfaces.  A  virtual  salt  flux  into  the  mixed  layer  will  occur  when  there  is  evaporation 
or  freezing  caused  by  decreasing  the  water  volume  while  keeping  the  salt  mass  relatively 
constant.  A  second  process  occurs  when  more  saline  water  from  below  is  mixed  into  the 
mixed  layer  by  entrainment. 

A  virtual  salt  flux  out  of  the  mixed  layer  occurs  when  there  is  precipitation  or 
melting  of  ice,  resulting  in  an  increase  in  water  volume  while  keeping  the  salt  mass 
relatively  constant.  Another  mechanism  occurs  when  there  is  entrainment  of  less  saline 
water  from  below  resulting  in  a  net  decrease  in  salinity.  Under  the  usual  conditions  in 
the  arctic,  the  water  under  the  mixed  layer  is  more  saline  so  this  condition  is  not  as 
common. 

The  salinity  of  the  water  affects  the  freezing  temperature  of  the  water.  The 
greater  the  salinity  the  lower  the  freezing  temperature  of  the  saline  water.  Hakkinen's 
(1987)  formula  for  the  freezing  temperature  of  salt  water  as  a  function  of  salinity  was: 

7)=  -  .003  -  .0527S  -  4.  x  10_5S2  (2.10) 

where 

•  Tf  is  the  freezing  temperature  (°  C) 

•  S  is  the  salinity  (parts  per  thousand) 

This  shows  a  nearly  linear  relationship  between  the  freezing  temperature  and  the  salinity 
due  to  the  S 2  term  being  small. 

5.  Salinity  and  Temperature  Effects  on  Buoyancy 

A  linearized  parameterization  for  density  at  the  water  surface  can  be  used  to 
demonstrate  the  relative  affects  of  salinity  and  temperature  fluxes  on  density: 

p  =  p0(l-jf(T-  T0)  +  J^(S~  S0))  (2.11) 

where: 

•  p  -  water  den  si  tv  (■—) 

m1 
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d  1 

— —  =  a  =  coefficient  of  thermal  expansion  of  water 
cl  L 

I 

-z~r  =  P  =  coefficient  of  salinitv  contraction  of  water  ( — — ) 
cS  PP1 

T  =  water  temperature  (°C) 

T0  =  reference  temperature  (°C) 

5  =  water  salinity  (ppt) 

Sa  =  reference  salinity  (ppt) 


The  parameters  a  and  p  are  function  of  temperature,  salinity,  and  pressure,  but  for  small 
changes  in  T  and  S,  a  and  /?  can  be  considered  constant. 

Equation  (2.11)  can  be  used  to  form  a  difference  equation  for  small  changes  to 
p,  giving: 


A  paz  —  a.pcAT  +  (ip0AS 


(2.12) 


g 

A  difference  buoyancy  equation  results  by  multiplying  equation  (2.12)  by  —  : 

r  O 

An 

Ab  =  g  =  -  cxgAT  +  figAS 

HO 


(2.13) 


where: 

•  Ab 

•  Ap 

•  8  = 

•  Po  : 

•  a.  = 

•  AT 

•  P  = 

•  AS 

Figure  5  shows  the  relative  values  of  ag  and  Pg  for  various  temperatures  and 
salinities.  In  the  extratropical  oceans  the  values  of  ag  and  Pg  are  on  the  order  of  25  x 
104  and  73  x  104  respectively.  When  a  slight  warming  of  the  surface  occurs,  Ab  is 
negative,  producing  a  stabilizing  decrease  in  surface  density.  In  the  polar  regions,  the 
temperatures  are  near  the  freezing  point  of  sea  water.  Figure  6  shows  an  expansion  of 


=  difference  in  buovancv 

kg 

=  difference  in  density  (~^f) 
=  acceleration  of  gravity  (-jj-) 
-  reference  density  (ppt) 


=  coefficient  of  thermal  expansion  (constant  for  small  changes)  (~Tq) 
=  difference  in  temperature  (°C) 

=  coefficient  of  salinity  contraction  ( — - — ^ 

=  difference  in  salinity  (ppt) 


PPt' 
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Figure  5  near  0°C  .  In  this  region  the  values  of  ag  and  (3g  are  on  the  order  of  2  x  10* 
and  78  x  10*  respectively.  A  slight  temperature  change  causes  a  negligible  change  in 
density,  but  a  slight  salinity  increase,  from  salt  rejection  during  the  freezing  of  salt 
water,  can  substantially  increase  the  density,  producing  a  destabilizing  increase  in 
density,  or  decrease  in  buoyancy  at  the  surface. 

6.  Entrainment 


The  mixed  layer  depth  will  increase  if  water  mass  comes  from  either  the  surface, 
ie.  rain,  snow,  or  melting  ice,  or  from  below.  When  the  mass  comes  from  below  by 
mixing,  it  is  called  entrainment.  Garwood  (1977)  and  others  have  defined  the 
entrainment  velocity  to  be 


IV. 


8h 

dt 


(2.14) 


where 

•  W,  -  the  entrainment  velocity 

•  h  -  the  mixed  layer  depth 

•  t  =  time 


B.  THE  TURBULENT  KINETIC  ENERGY  BUDGET 

1.  Total  Turbulent  Kinetic  Energy  Equation 

Garwood  (1987)  uses  the  total  turbulent  kinetic  energy  (TKE)  equation  to 
quantify  the  rates  associated  with  the  production,  transport,  storage,  and  dissipation  of 
turbulent  kinetic  energy.  This  equation  is  derived  from  the  Navier-Stokes  equation  for 
fluids  following  Garwood's  method.  The  total  TKE  equation  is 
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2.  Total  TKE  Change 


(2.15) 


The  term  on  the  left  side  of  equation  (2.15), 


describes  the  change  of  the  total  TKE  with  time.  The  forcing  functions  on  the  right  side 
of  the  equation  will  determine  how  the  total  TKE  changes  with  time.  This  term  can  be 
either  positive  or  negative,  for  an  increase  in  TKE  or  decrease  in  TKE  respectively. 

3.  Stress  Terms 


The  first  term  on  the  right, 


is  the  change  in  TKE  due  to  a  vertical  shear  in  the  horizontal  mean  flow.  This  term  is 
always  positive,  or  it  will  always  contribute  to  an  increase  in  turbulent  kinetic  energy. 
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When  a  wind  stress  in  the  +x  direction  acts  on  the  water  surface,  it  usually  will  set  up 
a  positive  gradient,  ~->0,  and  a  negative  Reynolds  stress,  u'w'  <0  ,  with  the  result 
that  an  increase  in  wind  stress  will  increase  the  TKE. 


4.  Buoyancy  Term 

The  second  term  on  the  right, 

w'p' 

~  ~S 

describes  the  change  in  TKE  due  to  a  vertical  density  flux,  or  equivalently  a  buoyancy 
flux.  This  term  will  decrease  TKE  (term  is  negative)  if  the  density  increases  with  depth, 
a  stable  condition.  An  unstable  condition  where  density  decreases  with  depth  would  be 
associated  with  an  increase  in  TKE  with  time.  Thermodynamic  forcing  or  salinity 
forcing  could  cause  an  unstable  situation;  such  conditions  could  be  heating  from  below, 
cooling  from  above,  evaporation  at  the  water  surface,  or  salt  rejection  from  the  freezing 
of  salt  water  at  the  surface.  All  these  situations  can  occur  in  freezing  regimes.  When 
this  term  is  negative  it  acts  to  dampen  TKE.  When  it  is  positive  it  acts  to  produce  TKE 
by  buoyant  production. 

5.  Turbulent  Transport  Terms 


The  third  term  on  the  right, 


consists  of  two  types  of  terms.  The  first  three  on  the  left 


redistribute  turbulence  in  the  vertical  while  the  forth  term 
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is  a  pressure  transport  of  TKE.  This  term  will  be  positive  or  negative  depending  upon 
whether  the  transport  causes  TKE  to  converge  or  diverge  from  a  given  level  in  the 
vertical. 


6.  Dissipation  Terms 
The  last  series  of  terms  on  the  right 
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are  the  dissipation  terms.  They  will  always  be  negative  and  will  tend  to  decrease  the 
TKE  with  time.  They  are  also  called  the  viscous  damping  terms,  meaning  the  turbulence 
is  dampened  out  due  to  the  viscosity  of  the  fluid.  This  is  the  principle  loss  term  for  the 
TKE  budget. 

7.  TKE  Summary 

The  Total  TKE  equation  explains  the  TKE  budget:  the  way  turbulent  kinetic 
energy  is  changed  with  time.  The  major  sources  of  TKE  are  the  horizontal  velocity  shear 
in  the  vertical  and  buoyancy  instabilities,  while  the  primary  sinks  of  TKE  are  increasing 
the  buoyancy  of  the  fluid  in  the  form  of  an  increase  of  potential  energy  and  the  viscous 
dissipation  off  TKE. 
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III.  THE  MIXED  LAYER  MODEL 


The  processes  involved  in  the  freezing  of  salt  water  are  not  completely  understood, 
and  there  are  various  ways  of  modeling  the  thermodynamics  of  the  freezing  process. 
Although  the  details  of  the  physics  are  not  completely  understood,  the  most  important 
processes  have  been  studied  extensively.  In  developing  a  model  that  is  be  used  to  study 
the  interaction  of  the  various  processes,  the  salient  processes  are  parameterized  using 
various  degrees  of  sophistication.  To  get  a  better  understanding  of  the  interaction  of  the 
relevant  variables,  it  is  best  to  start  with  a  simple  model  to  understand  fully  the  process 
interactions  before  developing  a  more  complex  model.  Although  some  systems  of  model 
equations  can  be  solved  analytically  and  valuable  information  can  be  gleaned,  such 
models  usually  have  been  simplified  and  may  lack  realism. 

For  this  research,  simplifying  assumptions  are  made  in  the  general  equations  to 
study  more  completely  the  various  parameters  of  mixed  layer  dynamics  and  their  effect 
on  the  onset  of  freezing.  To  evaluate  the  most  important  processes,  a  variation  of  the 
Kraus-Turner  one-dimensional  mixed  layer  model  (Kraus  and  Turner,  1967)  was  chosen 
and  simplifying  assumptions  were  made  to  investigate  the  relative  importance  of  the 
forcing  terms,  initial  conditions,  and  their  dynamic  interactions. 


A.  THE  ONE-DIMENSIONAL  MODEL 

The  model  used  for  this  study  is  derived  from  the  total  turbulent  kinetic  energy 
equation,  equation  (2.15),  with  simplifying  assumptions  and  simplifying 
parameterizations. 

I.  Model  Assumptions 


The  model  assumptions  following  Niiler  and  Kraus  (Niiler  and  Kraus,  1977)  are: 

1.  The  mean  temperature,  salinity,  and  horizontal  velocity  are  uniform  within  the 
mixed  layer. 

2.  On  the  depth  and  time  scales  of  the  model,  temperature  and  salinity  discontinuities 
can  exis.  across  the  lower  boundary  of  the  mixed  layer.  Thus  we  neglect  the  effect 
of  diffusion  and  conduction  across  this  lower  boundary. 
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3.  Temperature  changes  associated  with  the  frictional  dissipation  ofkinetic  energy  can 
be  neglected. 

4.  There  is  horizontal  homogeneity  of  wind  stress,  turbulent  kinetic  energy,  salinity, 
temperature,  and  viscous  dissipation.  Thus  horizontal  advection  is  neglected. 

5.  Only  the  case  of  a  deepening  mixed  layer  associated  with  a  freezing  regime  will  be 
considered.  Shallowing  mixed  layers  will  not  be  considered. 

6.  The  temperature  of  the  water  surface  is  approximately  the  temperature  of  the 
mixed  layer.  Thus  ice  formation  will  only  occur  when  the  mixed  layer  is  at  the 
freezing  point,  neglecting  the  effect  of  supercooling  at  the  water  surface. 

7.  The  effects  of  internal  waves  are  ignored. 


2.  Vertically  Integrated  TKE  Equation 

Using  the  above  assumptions  the  TKE  equation  is  integrated  vertically  over  the 
depth  of  the  mixed  layer.  First,  a  further  assumption  is  made  that  the  TKE  budget  is 
in  an  approximate  steady  state,  with  a  balance  between  production  and  dissipation. 
TKE  is  produced  by  both  buoyancy  flux  and  stress  production. 

The  first  term  of  equation  (2.1),  the  total  TKE  change  term  will  be  zero,  not 
allowing  the  total  TKE  to  change.  Thus, 


a.  Stress  Terms 

The  first  term  on  the  right  of  equation  (2.1),  the  shear  production  term,  is 
integrated  vertically,  giving 


where  is  the  friction  velocity.  The  friction  velocity,  u.  is  defined  as 
\w(uex  +  vey)\  =  —  m  u.  =  —  Caua 

where 

•  p  is  the  water  density 
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•  p '  is  the  air  density 

•  t„  is  the  wind  stress  magnitude  at  the  surface 

•  is  the  velocity  of  the  wind  at  a  reference  height  (normally  10  meters) 

•  C„  is  the  drag  coefficient  at  the  reference  level 

If  8u  is  proportional  to  u.,  then 


*o 


u  w 


cu 

8z 


a  ul 


(3.1) 


where  a  ~  6. 


b.  The  Buoyancy  Flux 

The  second  term  on  the  right  in  equation  (2.15)  can  be  rew'ritten  using  the 
definition  of  buoyancy, 


b 


Po~  P 
Po 


g 


Po-(p-  p') 
- To - * 


But, 


b 


Po~  P 
Po 


g 


and  thus 


b'  = 


P'g 

Po 


This  gives 


w  p 

~p7 


g  =  b'w' 


The  vertically  integrated  buoyancy  flux  thus  can  be  expressed  as 


b'w'dz 


h[b'w'(  0)  +  b’w'{-h)] 
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where 


•  h  is  the  mixed  layer  depth 

•  6'w'(0)  is  the  buoyancy  flux  at  the  surface 

•  b'w'(  —h)  is  the  buoyancy  flux  at  the  lower  interface 

Using  the  linearized  equation  of  state 

p  =  pji  1  -  <x(T-T0)  +  fi(S-Soy] 


where 

•  p  is  the  density 

•  p„  is  the  reference  density  at  T„  and  Sc 

•  a  is  the  coefficient  of  thermal  expansion 

•  P  is  the  coefficient  of  haline  contraction 

•  g  is  gravity 

•  T  is  the  water  temperature 

•  T„  is  a  reference  temperature 

•  S  is  the  water  salinity 

•  SQ  is  a  reference  salinity 

the  mean  density  is 

p  =  p0[l  -  a(T~T0 )  +  P(S-S0)] 
leaving  the  fluctuating  component 

p'  =  p0[  -  «r  +  PS'} 


This  gi-  es 


b'w'  =  - 


p  w 

Po 


g  =  agTw'  -  pgS'w' 


The  vertically  integrated  buoyancy  flux  can  now  be  written  as 
•o 


P 

J-h 


b'w'dz  =  h(txgTw'(  0)  —  PgS'w'(  0)) 


+  h{a.gT w'(  —h)  -  pgS'w'(-h) ) 


(3.2) 


(3.3) 
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that 


c.  Transport  Terms 

The  transport  terms  can  also  be  related  to  the  friction  velocity  cubed,  so 


"0 

•'-a 


(3.4) 


d.  Dissipation  Terms 

The  dissipation  term  has  nine  components: 


Kraus  and  Turner  (1967)  incorrectly  neglected  dissipation,  but  here  it  is 
assumed  to  be  of  the  form: 


(3.5) 


gives 


e.  Reduced  TKE  Equation 

Substituting  equations  (3.2),  (3.3),  (3.4),  and  (3.5),  into  the  TKE  equation 
0  =  a  u.3  +  bu3  +  cu.3  +  h(agTw'{0)  —  (igS'w'{  0)) 


+  h(agT  w'(  ~h)  -  PgS'w'{-h )) 


For  0, 
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0 


=  —(aw.  +  bw.  +  cw.)  +  <xgTw’(  0)  —  pgS'w'(  0) 


+  agTw'(-h)  —  pgS'w'(-h) 


The  three  constants  of  proportionality  (a,  b  and  c)  are  combined  into  one 
constant,  C„  so 

r  3 

I  i  i  i  Ci  w* 

—  (aw*  -I-  bw*  +  cw* )  =  — t — 
h  n 


This  gives  the  fmal  form  of  the  TKE  equation, 
0  - 


r  3 

C,w. 


+  (xgT iv'(0)  -  pgS'w'i 0) 


4-  agTw'(-h)  —  P  gS’w’(-h) 
Using  equation  (3.2),  an  equivalent  form  of  the  TKE  budget  is 


(3.6) 


C  u 

0  =  +  I/V(0)  +  iv(-/i) 


(3.6A) 


The  physical  meaning  of  the  terms  of  equation  (3.6)  are: 

1.  Turbulent  kinetic  energy  is  created  in  the  mixed  layer  by  the  wind  stress  at  the 
surface,  and  it  is  distributed  over  the  full  depth  of  the  mixed  layer 

2.  The  excess  TKE  from  the  wind  stress  is  moderated  by  the  buoyancy  flux  associated 
with  the  heat  flux  and  the  salinity  flux  at  the  surface. 

3.  The  remaining  TKE  is  dampened  by  the  negative  buoyancy  flux  at  the  bottom  of 
the  mixed  layer  by  entraining  denser  deep  water  from  below. 


3.  Buoyancy  Fluxes 


The  buoyancy  fluxes  at  the  surface  and  the  bottom  of  the  mixed  layer  will  be 
dependent  on  the  heat  fluxes  and  the  salinity  fluxes  at  the  top  and  bottom  of  the  mixed 
layer.  For  this  study  of  the  onset  of  freezing,  the  heat  flux  at  the  surface  will  be 
transformed  mostly  into  a  salinity  flux  via  the  thermodynamic  equations  (3.7,  3.15,  3.16, 
3.17,  and  3.18).  Depending  on  the  initial  conditions  below  the  mixed  layer,  there  may 
be  a  stability  problem,  resulting  in  the  overturning  of  the  whole  water  column.  These 
cases  will  not  be  investigated  here. 
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a.  Surface  Heat  Flux 

The  Tw^O)  term  represents  the  surface  heat  flux  into  the  mixed  layer, 


pCprw'( 0)  =  -  Qn 


(3.7) 


where 

•  Q„  is  the  net  downward  heat  flux  at  the  surface 

•  p  is  the  density  of  the  water 

•  Cp  is  the  specific  heat  of  water 

•  T'w'(0)  is  the  temperature  flux  associated  with  the  heat  flux 

In  a  freezing  situation,  which  is  the  emphasis  in  this  study,  there  exists  a 
state  when  any  further  heat  flux  out  of  the  mixed  layer  will  not  result  in  further  cooling 
but  will  result  in  the  formation  of  ice. 

b.  Surface  Salinity  Flux 

The  S'w'(O)  term  represents  the  surface  salinity  flux  out  of  the  mixed  layer, 
5V(0)  =  -  S(£  -  P)  -  (S  -  S;)(F-  M)  (3.8) 


where 

•  S  is  the  salinity  of  the  mixed  layer  in  g,'kg 

•  5,  is  the  salinity  of  the  ice  formed  or  melted  in  g  kg 

•  £  is  the  evaporation  rate  in  rn'sec 

•  £  is  the  precipitation  rate  in  m/sec 

•  £  is  the  freezing  rate  in  m/sec 

•  M  is  the  melting  rate  in  m/sec 

c.  Surface  Buoyancy  Flux 

The  surface  buoyancy  flux  can  be  wTitten  as  a  function  of  the  heat  flux  and 
the  components  of  the  salinity  flux,  using  equations  (3.2),  (3.7),  and  (3.8): 

*V(0)  =  +  ME  -  P)S  +  MF  -  ms  ■  S/)  (19) 


28 


d.  Heat  Flux  at  the  Lower  Boundary 

The  temperature  flux  at  the  bottom  of  the  layer  is  Tw'(  —h)  and  represents 
the  thermal  energy  brought  from  the  deeper  water  into  the  mixed  layer  by  entrainment. 
Here  the  effect  of  molecular  conductivity  of  thermal  energy  across  the  interface  is 
neglected.  This  allows  the  temperature  flux  to  be  written  as 

T V(-A)  =  —  AT  We  (3.10) 

where 

AT  =  T  ~  T0 


and  where 

•  h  is  the  mixed  layer  depth  in  meters 

•  T  is  the  mixed  layer  temperature  in  °C 

•  AT  is  the  temperature  jump  at  the  base  of  the  mixed  laver  between  the  mixed  layer 
in  °C 

•  T„  is  the  deep  water  temperature  in  °C 

•  \V,  is  the  entrainment  velocity  in  m/sec 

e.  Salinity  flux  at  the  Lower  Boundary 

The  salinity  flux  at  the  base  of  the  mixed  layer  S’w’(  —h)  ,  is  caused  by  the 
entrainment  of  lower  water  with  a  different  salinity.  This  can  flux  can  be  written  as 

SV(-A)  =  -  AS  We  (3.11) 

where 

AS  =  S  -  S0 


and  where 

•  h  is  the  mixed  layer  depth  in  meters 

•  AS  is  the  jump  salinity  between  the  mixed  layer  salinity  and  the  deep  water  salinity 
in  ppt 

•  S  is  the  mixed  layer  salinity  in  ppt 

•  S„  is  the  deep  water  salinity  in  ppt 

•  W,  is  the  entrainment  velocity  in  m/sec 
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f.  Buoyancy  Flux  at  the  Lower  Boundary 

The  buoyancy  flux  at  the  base  of  the  mixed  layer  can  be  written  from  the 
flux  equations  at  the  base  of  the  mixed  layer,  equations  (3.10)  and  (3.11),  and  equation 
(3.2) 

W{-h)  =  -Ab  H'e  (3.12) 


where 


Ab  =  agAT  —  fig  AS 


(3.13) 


4.  Entrainment  Velocity 

The  Entrainment  Velocity,  Wt,  can  be  obtained  from  equation  (3.6A),  where  a 
constant,  C2  ,  will  be  included  to  tune  this  simplified  model  to  fit  observations. 


0  = 


C,m* 


+  C2b'w’{  0)  +  b'w'i-h) 


Using  (3.13),  this  equation  becomes 


0  = 


C,m. 


+  C/>'  vv'(W  _  iv 

4.  S  •*  C 


or 


A  bWt  - 


J 

C,u. 


+  C2b'w'(  0) 


Solving  for  the  entrainment  velocity,  IV„ 


W  =  — — 

Ab 


C,M.  -  N 

+  C2b'w'(  0) 


or 
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(3.14) 


W. 


r-  3 

_ 1 _ C\u' 

agAT  —  fig  AS'  h 


+  - I - 

agAT  -  pgAS' 


+  Pg(E  —  P)S  +  Pg(F-  M)(S  -  Sd 


Furthermore,  it  is  required  that  W,  ^  0  ,  or  only  deepening. 

The  model  tuning  constants  C,  and  C2  are  used  to  allow  the  model  to  fit 
observations.  Chu  and  Garwood  (1988)  indicate  that  the  value  between  2  and  1  should 
be  used  for  C,,  and  a  value  near  .2  should  be  used  for  C,  .  For  this  study  a  normal  value 
of  2.  for  Ci  and  .2  for  C2  will  be  used. 


When  the  value  of  A b  nears  zero  the  model  becomes  computationally  unstable. 
This  condition  is  beyond  the  scope  of  this  study.  Thus  the  entrainment  velocity  here 
will  have  a  realistic  upper  limit,  consistent  with  observations. 


B.  THE  THERMODYNAMIC  MODEL 

The  thermodynamic  equations  will  relate  the  energy  needed  to  lower  the  water 
temperature  to  the  freezing  point  and  the  energy  removal  needed  to  freeze  the  water. 
Since  the  mixed  layer  regimes  for  the  onset  of  freezing  are  the  primary  focus  of  this 
study,  only  the  conditions  involving  a  net  upward  surface  heat  flux  are  of  interest. 

1.  The  Surface  Heat  Flux 

The  net  heat  fluxed  out  of  the  mixed  layer  is: 

Qn  =  Qr  +  Ql  +  Qe  +  Qs 


where 

•  =  net  heat  flux 

•  Qr  ~  heat  flux  from  shortwave  radiation 

•  Q,  =  heat  flux  from  longwave  radiation 


•  Q,  -  latent  heat  flux 

•  Q,  =  sensible  heat  flux 
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From  chapter  II,  equations  (2.1),  (2.5),  (2.6),  (2.7),  give 

Qn  =  Fc(l  -  kc3)(l  -a)  +  .7855(1  +  .2232c2'75)cr7^  -  .97*7? 
+  PaLCeU{qa  -  q0)  +  pacpCsL\Ta  -  Ts) 


where: 

•  Qn  =  net  heat  flux 

•  F„  =  incoming  solar  radiation 

•  k  =  seasonal  cloudiness  parameter 

•  c  —  cloudiness  in  tenths 

•  a  =  albedo  of  the  surface 

•  a  —  Stefan-Boltzman  constant 

•  T„  =  air  temperature  at  the  reference  level 

•  p,  =  air  density 

•  L  =  latent  heat  of  vaporization 

•  C,  =  latent  heat  bulk  transfer  coefficient 

•  U  -  wind  speed 

•  q0  -  specific  humidity  of  the  air  at  the  reference  level 

•  q„  -  specific  humidity  of  the  air  at  the  reference  level 

•  p„  -  density  of  the  air 

•  c„t  -  specific  heat  of  air 

•  C,  -  sensible  heat  bulk  transfer  coefficient 

•  T,  =  temperature  at  the  surface 

To  simplify  the  equations,  only  the  polar  night  will  be  investigated  thus 
eliminating  any  solar  radiation,  i.e.  Q,  =  0.  Other  assumptions  are  that  the  specific 

humidity  at  the  surface  will  be  at  saturation  and  the  temperature  at  the  surface  will  be 

the  mixed  layer  temperature,  Tr  —  T  .  The  values  of  C,  and  C,  are  assumed  constant. 
This  gives  an  equation  that  is  dependent  on  the  air  temperature,  specific  humidity  of  the 
air,  sea  surface  temperature,  wind  speed  and  cloudiness: 

Qn  =  .7855(1  +  . 2232c27V^  +  paLCeU{qa- q0) 

+  pacpCsU(Ta-T)  -  .91  oT* 
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Figure  7.  Heavyside  Step  Function  And  Modified  Heavyside  Step  Function:  The 
Heavyside  Step  Function  (solid  line)  which  is  0  for  values  of  T—  Tf<,  0 
and  is  1  for  values  of  T  —  Tf  >  0  .  This  function  has  sharp  corners  and 
is  second-order  non-continuous  at  0.  The  dotted  line  is  a  modified 
Heavyside  Step  Function  which  has  smooth  comers  and  is  second-order 
continuous.  This  function  was  used  in  place  of  the  Heavyside  Step 
Function  to  reduce  computational  noise. 
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2.  The  Freezing  Temperature 


The  freezing  temperature  formula  of  Hakkinen  (1987)  is  used: 

7}  =  -.003  -.05275  -4.0xl0_552  (3.16) 


where 

•  Tf  is  the  freezing  temperature  of  salt  water 

•  5  is  the  salinity  of  the  water  (°%) 

To  ensure  no  loss  of  energy  when  the  mixed  layer  is  near  the  freezing 
temperature,  a  distinction  was  made  between  a  decrease  in  water  temperature  and  the 
production  of  ice  using  the  Heavyside  Step  function,  A: 


A (r,  7})  = 


0  for  T  <  Tf 
1  for  T>  Tf 


Figure  7  shows  the  Heavyside  Step  Function. 


To  reduce  computational  noise  caused  by  this  second  order  non-continuous 
function,  A  is  approximated  with: 

1  +  TANH((T  -  77)2200) 

MT,  Tf)  =  - - - L -  (3.17) 

Figure  7  shows  the  modified  Heavyside  Step  Function.  Here  the  corners  are 
slightly  rounded,  and  the  function  is  second  order  continuous.  This  function  eliminated 
computational  noise. 


3.  The  Freezing  Rate 


The  freezing  rate  was  determined  by  the  amount  of  energy  remaining  after  first 
lowering  the  mixed  layer  temperature  to  freezing: 

F  =  ~tf{~  ~h;  +  ATW^[  -  A(r-7>))  (3-18) 


34 


where: 


•  F  is  the  freezing  rate 

•  cp  is  the  specific  heat  of  water 

•  Lf  is  the  latent  heat  of  fusion  of  ice 

•  Qn  is  the  net  heat  flux 

•  p  is  the  water  density 

•  AT  is  the  difference  between  the  mixed  layer  and  lower  layer  temperature 

•  Wt  is  the  entrainment  velocity 

•  A  is  the  Heavyside-step  function  approximation 

C.  THE  NUMERICAL  MODEL 

In  the  numerical  model,  ordinary  differential  equations  are  used  of  the  form: 

4r  =  {[-iivw  - 1  -  nv(-€  <3-19) 

1.  The  Time  Dependent  Equations 

From  equation  (3.19)  the  time  dependent  thermal  equation  is 
=  j[-Tw(0)  -  {  -  7V(  -h)}]A(T,  Tf) 

Using  equations  (3.7)  and  (3.10)  this  can  be  written  as 

17  -  i  %;  -ATW.MT.Tj>  (3.20) 

Similarly  the  time  dependent  Salinity  equation  is 

~  =  I[-SV(0)-{-SV(-A)}] 

Using  equations  (3.8)  and  (3.11)  this  can  be  written  as 

=  j[S(£-P)  +  (S-SJiF-M)  -  ASWj  (3.21) 
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The  time  dependent  equation  for  the  buoyancy  is 


-f-  -  j-[-f>V(0)  -  {  -  4V(-A)}] 

From  equations  (3.9)  and  3.12)  this  can  be  written  as 

8b  _  1 

dt  h 

The  time  dependent  mixed  layer  depth  equation  is 

~  =  We-w{-h)  -  (E  -  P  +  F  -  M) 


ag 


A 

PCp 


-  pg(E  -  P)S  -  pg(F-  mS  -  S,)  -  A b  We 


(3.22) 


(3.23) 


where 

•  IV.  is  the  entrainment  velocity 

•  h(  —  h)  is  a  time  dependent  function  for  vertical  velocity  (i.e.  internal  wave,  tides) 

•  £  is  the  evaporation  rate 

•  £  is  the  precipitation  rate 

•  £  is  the  freezing  rate 

•  M  is  the  melting  rate 


The  time  dependent  equation  for  ice  growth  is 


(3.24) 


where 

•  h,  is  the  thickness  of  the  ice 

•  £  is  the  freezing  rate 

•  p  is  the  density  of  the  w'ater 

•  p,  is  the  density  of  the  ice 

2.  Model  Formulas 

All  the  formulas  used  to  define  the  model  variables,  and  the  model  constants, 
are  summarized  below  and  in  Tables  I  through  4. 


36 


1  clu 

lVg  =  agAT  -pgAS'  h 


(3.14) 


agAT  -  fig  AS 


r  Qn 

~  C 2  ~  aS  pcp 


+  pg(E-  P)S  +  pg(F-  M){S 


-Sd 


for  IV,  ;>  0,  otherwise  W,  =  0 


‘p  (  Qn 


F~  Lf[  Pcp 


+  ATlVep  -  A (T,Tf)) 


(3.18) 


b'W'(  0)  =  +  Ps(E~P)s  +  Pg(F-M){S-Si) 


b'w'i-h )  =  -A b  lVe 


(3.12) 


Ab  =  agAT  —  fig  AS 


(3.13) 


AT  =  T  -  Tn 


AS  =  S  -  Sn 


«•  = 


Tf  =  -.003  — .0527S  -  4.0xl0_s52 


(3.16) 


1  +  TANH((F  -  77)2200) 
A(T,  Tf)  -  - 2 - - 


(3.17) 


0*  -  .7855(1  +  ,2232c2'75)ct7^  +  paLCeU{qa-  q0) 
+  pac  CsU(Ttt-T)  -  Slot 


(3.15) 
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Table  1.  MODEL  VARIABLES 


Variable 

Units 

Description 

a 

1 

°C 

Coefficient  of  Thermal  Expansion 

P 

kg 

g 

Coefficient  of  Salinity  Contraction 

b 

m 

s» 

Buoyancy  of  the  Mixed  Layer 

A  b 

m 

s2 

Buoyancy  Difference  between  the  Mixed  Layer  and  the  Deep 
Layer 

b’w'{ 0) 

m2 

s3 

Buoyancy  Flux  at  the  Surface  of  the  Mixed  Layer 

b'w'{  -h) 

m2 

s3 

Buoyancy  Flux  at  the  Bottom  of  the  Mixed  Layer 

c 

Tenths 

Cloudiness 

E 

m 

s 

Evaporation  Rate 

F 

m 

s 

Freezing  Rate 

h 

m 

Depth  of  the  Mixed  Layer 

h, 

m 

Thickness  of  the  Ice 

A 

Heavyside  Step  Function  (Modified) 

M 

m 

s 

Melting  Rate  of  the  Ice 

P 

m 

s 

Precipitation  Rate  into  the  Mixed  Layer 

<h 

Specific  Humidity  of  the  Air  at  the  reference  level 
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Table  2.  MODEL  VARIABLES  CONTINUED 


Variable 

Units 

Description 

q, 

Specific  Humidity  of  the  Air  at  the  Surface 

a 

W 

m2 

Net  Heat  Flux  downward  at  the  Surface  of  the  Mixed  Layer 
(Downward  is  Positive) 

p 

msm 

Density  of  the  Mixed  Layer  Water 

Pc 

3  1  ^ 

Density  of  the  Air  at  the  reference  level 

p, 

kg 

m3 

Density  of  the  Ice 

s 

g 

ks 

Salinity  of  the  Mixed  Layer 

S, 

g 

ke 

Salinity  of  the  Ice 

Sc 

g 

ke 

A  Salinity  of  the  Deep  Water 

AS 

g 

ke 

Salinity  Difference  between  the  Mixed  Layer  and  the  Deep 
Water 

t 

s 

Time 

T 

Temperature  of  the  Mixed  Layer 

Tc 

Temperature  of  the  Air  at  the  reference  level 

Tf 

Temperature  of  Freezing 

L 

Temperature  of  the  Deep  Water 

T, 

Temperature  of  the  Air  at  the  Surface  of  the  Water  (assumed 
to  be  equal  to  T) 
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Table  3.  MODEL  VARIABLES  CONTINUED 


Variable 

Units 

Description 

AT 

D 

Temperature  Difference  between  the  Mixed  Layer  and  the 

Deep  Water 

u. 

m 

s 

Friction  Velocity  of  the  Mixed  Layer 

U ,  u 

m 

s 

Velocity  of  the  Air  at  the  reference  level 

w 

m 

s 

A  Vertical  Velocity  Function  at  the  Interface  of  the  Mixed 

Layer  and  the  Deep  Water  due  to  some  Forcing  (i.e.  tides, 
internal  waves,  etc.) 

W. 

m 

s 

Entrainment  Velocity 
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Table  4.  MODEL  CONSTANTS 


Constant 

Units 

Value 

Description 

CP 

J 

kg°C 

4.18jc103 

Specific  Heat  of  Water 

J 

kg°C 

1004 

Specific  Heat  of  Air 

c. 

none 

1.3jdO-3 

Drag  Coefficient  for  the  Air  at  a  Reference 

Level 

c. 

none 

1.28jc10-3 

Latent  Heat  Bulk  Transfer  Coefficient 

C, 

none 

1.4jc10  3 

Sensible  Heat  Bulk  Transfer  Coefficient 

Q 

none 

2.0 

Tuning  Constant  for  Friction  Velocity  Term 

Q 

none 

0.2 

Tuning  Constant  for  Surface  Buoyancy  Flux 

Term 

g 

m 

s2 

9.83 

Acceleration  due  to  Gravity 

L 

J 

kg 

2.5jrl06 

Latent  Heat  of  Vaporization 

Lf 

J 

kg 

302.0jc103 

Latent  Heat  of  Fusion  of  Ice 

a 

kg 

s3(°K)J 

5.67x10-* 

Stefan-Boltzman  Constant 
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IV.  MODEL  RESULTS 


A.  MODEL  INITIAL  CONDITIONS  AND  BOUNDARY  CONDITIONS 

For  this  study  the  model  was  run  for  various  initial  conditions  and  boundary 
conditions.  For  boundary  conditions,  the  wind  speed  and  various  parameters  affecting 
Qn,  i.e.  air  temperature,  specific  humidity  of  the  air,  and  cloudiness,  were  considered. 
For  initial  conditions,  or  pre-conditioning  situations,  consideration  was  given  to  the 
initial  mixed  layer  salinity,  which  determines  the  mixed  layer  freezing  temperature,  the 
deep  water's  temperature,  the  deep  water's  salinity,  and  the  initial  mixed  layer  depth. 


Table  5.  NORMAL  VALUES  AND  SENSITIVITY  RUN  VALUES 


Variable 

Normal 

Values 

Other  Values  Used  For  Sensitivity  Analysis 

Wind  Speed, 

U 

5  — 

0,  .5,  1,  1.5.  2.  2.5,  3,  3.5,  4,  4.5,  5,  5.5.  6,  6.5,  7,  7.5,  8,  8.5, 

9,  9.5,  10,  10.5,  11,  11.5,  12 

Salinity 
Jump,  AS 

5  _£_ 

‘  kg 

.3,  .4,  .5,  .6,  .7,  .8,  .9,  1.0,  1.2,  1.5,  2.0,  2.5,  3.0,  4x,  5.0,  6.0, 
6.5,  7.0,  7.5,  8.0,  8.5,  9.0,  9.5.  10.0,  11.0,  12.0,  13.0,  14.0, 

15.0,  16.0.  17.0,  18.0,  19.0,  20.0,  25.0 

Temperature 
Jump,  AT 

2. 8688 
°C 

-1,  -.5,  0,  .5,  1.  1.5,  2,  2.2,  2.4,  2.6,  2.7,  2.8,  2.868,  2.9,  3.0, 

3.5,  4.0,  5.0,  10.0 

Mixed  Layer 
Depth,  h 

60  m 

5,  10,  30,  40,  50.  60,  70,  80,  90,  100,  120,  150,  160,  175,  200, 
250,  300,  500,  1000 

Surface  Heat 
flux,  Qn 

W 

-350  — 7 
m2 

100,  0, -100, -200, -300, -350, -400,  -450, -500,  -600,  -700, 

-800 

Model  simulations  were  made  with  the  mixed  layer  temperature  close  to  its  freezing 
temperature,  allowing  the  mixed  layer  to  begin  freezing.  The  specific  ranges  of  boundary 
conditions  ranged  from  the  most  severe  arctic  storm  conditions,  wind  >  60  knots,  to 
calm  winds,  while  varying  the  factors  affecting  Qn  between  physically  realistic  maximum 
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Wind  Velocity  (U)  ( ■y- ) 


Figure  8.  Freezing  Rate  as  a  Function  of  Wind  Speed,  U,  Only:  Results  of  model 
runs  for  variations  in  wind  speed,  without  changing  any  other  initial 
condition  or  boundary  condition.  For  a  constant  net  heat  flux  upward, 
there  is  a  maximum  freezing  rate  under  no*wind  conditions  where  the 
effects  of  wind  induced  turbulent  mixing  is  negligible.  Under  higher  wind 
conditions,  wind  speeds  over  1 1  the  heat  flux  into  the  mixed  layer 
from  the  entraining  of  warmer  water  belowr  overcomes  the  heat  flux 
upward  out  of  the  mixed  layer  and  net  mixed  layer  warming  occurs. 
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Figure  9.  Freezing  as  a  Function  of  u  *:  Results  of  model  runs  where  the  wind 
speed  was  varied,  without  varying  other  initial  conditions  or  boundary 
conditions.  The  results  show  the  linear  relationship  between  the  freezing 
rate  and  W„  which  at  small  freezing  rates,  is  a  linear  function  of  uJ. 
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and  minimum  values.  The  initial  salinity  and  water  temperature  values  were  varied 
between  realistic  limits.  Unrealistic  conditions  were  also  considered  to  determine  if  the 
model  would  produce  unstable  conditions. 

A  major  consideration  for  this  thesis  is  the  regimes  where  freezing  occurs  and  does 
not  occur,  and  the  boundary  between  them.  For  this  reason  the  freezing  rate  is 
considered  as  the  model  output  that  will  be  used  to  compare  the  results  of  the  model 
runs. 

Table  5  shows  the  values  that  were  considered  as  standard  initial  conditions  and 
boundary  values,  plus  those  values  that  were  used  on  the  model  runs. 

1.  Model  Sensitivity  To  Wind 

The  model  was  run  under  varying  ranges  of  wind  conditions  without  varying  any 

other  factors.  This  takes  into  consideration  regimes  where  the  wind  speed,  U,  increases 

or  decreases  with  an  appropriate  change  in  sensible  and  latent  heat  fluxes  to  keep  the 

IV 

net  heat  flux,  Q„,  constant.  As  an  example,  for  a  constant  condition  of  Q„  =  —350  — -, 

PH  ^ 

and  U  =  5—  ,  an  appropriate  air  temperature  would  be  -21  °C.  T„  —  —  21°C  .  For 
a  wind  speed  of  U  =  30 -y-  an  appropriate  air  temperature  would  be  -.03  °C. 
T,  =  -  ,03°C 

The  only  physically  unrealistic  situation  was  for  calm  winds,  where  the 

W 

maximum  attainable  realistic  Q.  is  on  the  order  of  -200  — r. 

ml 

Figure  8  shows  the  results  of  the  model  runs.  Because  only  the  wind  speed  was 

varied,  this  mostly  represents  the  effect  of  varying  the  entrainment  of  warmer  water  from 

below.  From  equation  (3.14),  if  Qn  is  constant  and  F  is  small,  the  entrainment  velocity 

acts  as  a  cubic  function  of  From  equation  (3.18),  again  if  Qn  is  constant,  the  freezing 

rate  should  be  a  function  of  &TW„  which  is  negative  under  these  conditions,  and  the 

resulting  curve  should  reflect  a  negative  turning  cubic  function.  Under  conditions  of  Q„ 

IV 

being  held  constant  at  -350  — r,  the  heat  flux  at  the  bottom  of  the  mixed  laver,  from 
entraining  warmer  water,  overcomes  the  heat  flux  out  of  the  surface  when  the  wind 
speed  nears  1 1  For  higher  wind  speeds  a  net  mixed  layer  warming  occurs.  Figure 
9  shows  the  relationship  of  the  freezing  rate  as  a  function  of  ul,  vice  U.  This  shows  the 
near  linear  relationship  of  the  W„  a  function  of  u.\  and  the  freezing  rate. 
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Figure  10.  Freezing  Rate  as  a  Function  of  the  Salinity  Jump:  Results  of  model 
runs  where  the  salinity  jump,  AS,  was  varied,  without  varying  other 
initial  conditions  or  boundary  conditions.  The  results  show  a  rapid 
decrease  in  freezing  rate  as  the  salinity  jump  nears  -.3,  and  a  slow  near 
linear  increase  for  AS  <  —2.5, 
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Figure  11.  Freezing  Rate  as  a  Function  of  Salinity  Jump:  An  extension  of  Figure 
10  for  values  of  Salinity  Jump,  AS,  to  -25.0  °C.  The  extended  values 
past  -5.0  show  the  leveling  off  trend  for  greater  magnitudes  of  AS  and 
its  approach  toward  the  maximum  freezing  available  of 
F  =  112.8x10-'. 
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Figure  12.  Freezing  Rate  as  a  Function  of  Temperature  Jump:  Results  of  model 
runs  where  the  temperature  jump,  AT  is  varied,  without  varying  the 
other  initial  conditions  or  boundary  conditions.  The  results  show'  a 
near  linear  decrease  in  freezing  rate,  F,  as  the  temperature  jump 
increases  in  magnitude.  For  values  of  AT  >  —1.8  ,  the  freezing  rate 
is  greater  than  would  be  expected  form  just  the  surface  thermal  flux 
alone.  As  the  magnitude  of  AT  nears  -10.0,  the  thermal  flux  from 
entrainment  is  greater  than  the  surface  thermal  flux  and  net  wanning 
of  the  mixed  layer  occurs  and  freezing  stops. 


2.  Model  Sensitivity  To  Salinity  Jump 


Figure  10  shows  the  results  of  only  varying  the  salinity  of  the  deep  layer, 
resulting  in  a  change  in  the  salinity  jump,  AS,  between  the  layers.  The  general  trend  of 
the  results  are  an  increasing  freezing  rate  with  an  increasing  salinity  jump,  or  increasingly 
negative  AT. 

For  conditions  where  the  salinity  of  the  two  layers  is  too  close,  AT  >  —.3,  there 
was  an  over  turning  of  the  water  column  and  no  freezing  occurred.  This  occurs  when 
the  mixed  layer  density  is  increased  by  increasing  the  salinity  to  a  point  where  its  density 
is  the  same  or  greater  than  the  deep  layer. 

When  the  salinity  jump  is  increased,  the  freezing  rate  begins  to  level  off.  This 
occurs  due  to  the  decreasing  heat  flux  at  the  base  of  the  mixed  layer.  From  equation 
(3.10),  the  thermal  flux  at  the  base  of  the  mixed  layer  is  a  function  of  the  entrainment 
velocity,  W„  and  the  temperature  jump,  AT .  The  temperature  jump  is  not  changing  but 
the  entrainment  velocity  is.  From  equation  (3.14),  the  entrainment  velocity,  W„  is  a 
function  of  the  inverse  of  the  buoyancy  difference,  A b.  As  the  salinity  difference 
increases  the  buoyancy  difference  increases  and  the  entrainment  velocity  decreases. 
From  equation  (3.18),  when  the  entrainment  velocity  decreases,  the  freezing  rate 
increases. 

Figure  1 1  shows  the  results  of  increasing  the  density  jump  to  much  higher  levels. 
This  shows  the  leveling  off  of  the  freezing  rate  at  higher  density  jump  values.  The  limit 
for  this  is  the  maximum  freezing  available,  near  112.8xlO-8  -y-,  the  situation  where 
there  is  relatively  no  entrainment  of  warmer  water  and  the  full  surface  heat  flux  is 
converted  into  the  production  of  ice.  The  large  density  difference  between  layers  acts  to 
virtually  insulate  the  two  layers  from  each  other.  This  relative  salinity  difference  can 
occur  when  fresh  water  runoff  overlies  more  saline  deeper  water  which  has  not  been  well 
mixed  by  wind  stress. 

3.  Model  Sensitivity  To  Temperature  Jump 

Figure  1 2  shows  the  results  of  varying  only  the  temperature  jump  by  varying 
only  the  deep  water  temperature.  The  general  trend  is  for  an  increase  in  freezing  rate 
with  a  lowering  of  the  deep  layer  temperature,  or  AT  gets  more  positive.  The  curvature 
of  this  curve  is  much  more  linear  than  the  curves  in  Figures  10  and  1 1  because  the  effect 
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Figure  13.  Freezing  Rate  as  a  Function  of  Mixed  Layer  Depth:  Results  of  model 
runs  where  the  mixed  layer  depth,  h,  is  varied,  without  varying  the  other 
initial  conditions  or  boundary  conditions.  The  results  show  a  rapid 
increase  in  freezing  rate,  F,  as  depth  increases  past  10  m  to  150  m.  This 
is  from  the  decreasing  effect  of  entrainment  of  warmer  deep  water  from 
wind  generated  TKE  as  the  depth  increases. 


0-0  2.0  4.0  6.0  B.O  10.0 


Figure  14.  Freezing  Rate  as  a  Function  of  Surface  Heat  Flux:  Results  of  model 
runs  where  the  surface  heat  flux,  Q„  is  varied,  without  varying  the  other 
initial  conditions  or  boundary  conditions.  The  results  show  the  near 


linear  increase  in  freezing  rate  with  an  increase  in  upward  heat  flux  out 
of  the  mixed  layer. 
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of  the  temperature  differences  on  Ab  is  so  much  less  than  the  effect  of  a  salinity 

79  0 

difference  from  the  small  value  of  a.  relative  to  B  of  i  ‘ 

2. 36 

In  the  region  of  AT  >  0,  the  freezing  rate  is  greater  than  the  maximum 
attainable  from  just  the  surface  heat  flux.  This  could  occur  when  a  more  saline  deep 
layer  is  near  its  freezing  temperature,  which  is  colder  than  the  mixed  layer  freezing 
temperature,  and  any  entrainment  of  deep  water  actually  increases  the  freezing  rate. 
This  can  be  seen  from  equation  (3.18),  when  AT  >  0,  the  two  thermal  fluxes  work 
together  to  increase  the  freezing  rate. 

At  the  other  extreme,  there  is  a  AT  so  great  that  entrainment  of  the  warmer 
water  overpowers  the  surface  heat  flux  and  a  warming  of  the  mixed  layer  occurs.  From 
equation  (3.18)  and  under  these  model  conditions  of  AT  <  —  10.0  ,  the  effect  of 
AT  W„  which  is  negative,  is  greater  in  magnitude  than  the  surface  thermal  flux,  —S~  , 
which  is  positive. 


4.  Model  Sensitivity  To  Mixed  Layer  Depth 

Figure  13  shows  the  results  on  the  freezing  rate  of  only  varying  the  mixed  layer 
depth.  The  general  trend  is  for  a  rapidly  increasing  freezing  rate  when  increasing  the 
mixed  layer  depth  from  3  m  to  50  m,  a  gradual  increase  of  freezing  rate  for  mixed  layer 
depths  from  50  m  to  175  m  ,  and  then  a  near  asymptotic  increase  in  the  freezing  rate  for 
mixed  layer  depths  greater  than  175  m. 

The  rapidly  increasing  freezing  rate  for  small  values  of  mixed  layer  depth  reflect 

the  relative  importance  of  the  TKE  generated  by  wind  stress  distributed  over  the  mixed 

1  CTj  w? 

layer.  From  equation  (3.14),  when  h  is  small,  the  — - 7—  term  is  very  large  and  varies 

Ab  n 

inversely  with  mixed  layer  depth,  h.  This  results  in  a  large  value  for  the  entrainment 
velocity,  IV„  which  rapidly  overpowers  the  thermal  flux  at  the  surface. 

j  C  w? 

As  the  mixed  laver  depth  increases  the  relative  value  of  the  — - r~  term  is 

Ab  h 

increasingly  smaller  and  becomes  less  important  in  the  freezing  of  the  mixed  layer  and 
the  freezing  rate  approaches  that  of  the  maximum  possible  freezing  rate.  This  is 
reflected  in  the  figure  by  the  assymtotic  approach  to  a  constant  value  as  h  increases. 


5.  Model  Sensitivity  To  Surface  Heat  Flux 


Figure  14  shows  the  nearly  linear  relationship  of  the  surface  heat  flux,  Q„  and 
the  freezing  rate.  The  general  trend  is  a  linear  increase  in  freezing  rate  as  the  surface 
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heat  flux  increases.  When  the  surface  heat  flux  is  close  to  0  or  positive  there  is 
insufficient  cooling  to  counter  the  entrainment  flux  that  causes  a  net  warming  of  the 
mixed  layer. 

From  equation  (3.14),  the  entrainment  velocity  is  a  nearly  linear  function  of  Qn 
since  h  and  u.3  are  held  constant.  From  equation  (3.18),  the  freezing  rate  should 
approximate  a  linear  function  since  IV,  will  be  a  nearly  linear  function  of  Qn,  and  AT  is 
held  constant. 


B.  NON-DIMENSIONAL  PARAMETERIZATION 


1.  Formation  of  Non-Dimensional  Parameters 


To  understand  the  relationship  between  many  varying  variables  it  is  useful  to 
organize  the  variables  into  non-dimensional  parameters  that  can  give  significant 
information  when  analyzed.  Non-dimensional  parameters  can  be  formed  from  the 
various  forcing  variables  and  boundary  variables  that  will  simplify  the  understanding  of 
the  mixed  layer  interaction  problem. 

2.  Freezing  Efficiency  Parameter 


The  first  parameter  is  a  freezing  efficiency  parameter.  This  is  the  ratio  of  the 

freezing  rate  to  the  heat  loss  from  the  surface,  or  the  maximum  heat  loss  normally 

available  for  freezing.  The  parameter  will  be  called  £*,  for  Efficiency.  To  formulate  the 

E  parameter,  the  freezing  rate,  F,  in  — ,  will  be  divided  bv  the  maximum  freezing 
Qn  m 

available,  ,  which  is  also  in  — .  Therefore  E  is  defined  as 
Lffi 


* 

£  = 


F 


Oil 

LfP 


To  make  the  E  positive  under  normal  conditions,  the  negative  of  Qn  will  be  used. 


E  = 


FPLf 

-Qn 


(4.1) 
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This  parameter  has  a  value  of  unity  when  all  the  energy  given  up  as  heat  flux  from  the 
surface  is  used  in  the  production  of  ice.  Whenever  this  parameter  is  less  than  1.0,  some 
of  the  heat  flux  from  the  surface  has  been  countered  by  a  heat  flux  into  the  mixed  layer 
by  entrainment  of  warmer  water  from  below.  The  only  way  the  absolute  value  of  this 
parameter  should  be  greater  than  1.0,  would  be  when  a  less  saline  mixed  layer  overlies 
a  more  saline  deeper  layer  that  is  near  freezing  temperature,  and  entrainment  cools  the 
mixed  layer.  This  causes  a  negative  heat  flux  across  the  lower  boundary  of  the  mixed 
layer,  and  it  would  give  a  value  of  E  >  1.0  . 

3.  Forcing  Ratio  Parameter 


A  ratio  of  the  wind  forcing,  which  provides  TKE  for  mixing  and  entrainment, 
to  the  buoyancy  flux  at  the  surface,  also  providing  TKE  for  mixing  and  entrainment, 
would  provide  information  about  the  relative  magnitudes  of  the  two  TKE-generating 
factors.  This  parameter  will  be  called  F,  for  forcing,  and  can  be  defined  as: 


or 


«.3 


h 

F^(0) 


hb'w'(0) 


(4.2) 


When  F  <  1.0  the  buoyancy  flux  from  freezing  dominates  the  wind  stress  at  the 
surface.  This  can  occur  when  the  mixed  layer  is  deep,  the  wind  is  low,  or  rapid  freezing 
occurs  from  a  large  heat  flux  to  the  atmosphere. 

When  F  >  1.0  then  the  wind  forcing  dominates  the  buoyancy  flux  at  the  surface. 
This  occurs  when  the  wind  speed  is  high,  the  mixed  layer  is  shallow,  or  there  is  little 
freezing  due  to  entrainment  of  much  warmer  water  or  due  to  very  little  heat  flux  from 
the  mixed  layer  to  the  atmosphere. 
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4.  Stability  Ratio  Parameter 


A  ratio  of  the  two  parameters  which  affect  the  stability  of  the  mixed  layer  over 
the  deeper  layer  give  information  about  the  overall  stability  of  the  two  systems.  Under 
normal  freezing  conditions  in  the  Arctic,  a  cooler  layer  over  a  warmer  layer  decreases 
stability  or  buoyancy  of  the  mixed  layer.  A  less  saline  layer  over  a  more  saline  layer 
increases  stability  or  buoyancy  of  the  mixed  layer.  The  ratio  is  a  measure  of  the 
effectiveness  of  the  salinity  jump  in  maintaining  stability  over  the  temperature  jump's 
effectiveness  in  decreasing  stability.  The  stability  ratio  will  be  called  S’,  and  can  be 
formed  by  dividing  the  buoyancy  due  to  the  temperature  jump  over  the  buoyancy  due 
to  the  salinity  jump: 


pgAS  pAS 


(4.3) 


When  S'  =  1.0  ,  the  two  layer  are  neutrally  buoyant,  meaning  there  is  no 
density  difference  between  the  two  layers  even  though  there  may  be  a  salinity  and 
temperature  difference  between  the  layers.  In  this  case,  the  model  is  computationally 
unstable  and  would  predict  the  complete  overturning  of  the  wTater  column. 

When  S'  >  1.0,  the  temperature  jump  is  dominant  and  entrainment  may  lead 
to  a  warming  of  the  mixed  layer  and  may  stop  freezing  altogether. 

A  density  ratio  has  been  used  by  Schmitt  (1918,  1987,  1988)  and  others  to 
described  domains  where  thermohaline  staircases  could  be  found.  Schmitt  (1988)  defines 
his  ratio  the  density  ratio,  Rp  ,  but  used  T2  -  vertical  temperature  gradient,  instead  of 
the  temperature  jump,  AT,  and  S,  =  vertical  salt  gradient,  instead  of  the  salinity  jump, 
AS.  From  his  research  many  examples  of  thermohaline  staircases  have  been  found  and 
documented.  He  has  found  that  staircases  are  most  likely  when  1.0  >  R„  >  1.6  . 

When  S’  1.0  ,  the  salinity  jump  dominates  and  entrainment  would 
contribute  very  little  to  a  heat  flux  from  the  deep  layer  to  the  mixed  layer.  The 
entrainment  salinity  flux  into  the  mixed  layer  would  have  but  a  slight  impact  by  lowering 
the  freezing  temperature.  This  condition  would  lead  to  the  most  efficient  freezing 
situations. 


55 


Table  6.  NON-DIMENSIONAL  PARAMETERS 


Non-Dime 

Variable 

nsional  Formula 

Significance 

£* 

pLfF 

-Q. 

F 

Q, 

L  pLf  J 

Efficiency  of  freezing.  Actual 
freezing  rate  over  the  thermal  energy 
lost  from  the  mixed  layer  that  is 
available  for  freezing.  A  value  near 
one  indicates  very  little  energy  was 
entrained  from  the  deep  layer. 

F 

ui 

hb'w\ 0) 

-  ul  - 

h 

b'w\ 0) 

Turbulence  Forcing  Ratio. 
Turbulence  forcing  from  wind  stress 
over  turbulence  forcing  from  surface 
buoyancy  flux. 

S' 

a  AT  agAT 

PAS  PgAS 

Stability  ratio.  Ratio  of  buoyancy 
change  due  to  a  change  in 
temperature  over  the  buoyancy 
change  due  to  a  change  in  salinity. 

The  smaller  the  value  of  the 
parameter  the  greater  the  stability 
between  lavers. 

When  S‘  <  0  ,  the  temperature  and  salinity  density  differences  both  contribute 
to  decreasing  the  buoyancy.  The  only  physically  possible  case  would  be  to  have  warmer 
and  less  saline  water  over  cooler  and  more  saline  water.  If  the  mixed  layer  were  near 
freezing,  any  entrainment  could  increase  the  freezing.  This  is  the  situation  mentioned 
above  for  the  case  when  F  <  0.0  also.  This  could  occur  when  fresh  water  run  off 
encounters  calm  very  saline  salt  water  close  to  the  coast  and  spreads  out  as  a  thin  layer 
over  the  salt  water.  The  physical  process  of  conduction  has  been  ignored  in  this  model, 
which  would  be  the  major  heat  transfer  process  in  this  situation,  and  this  situation  would 
not  be  handled  well  by  this  model.  Table  6  summarizes  the  non-dimensional  parameters 
considered  in  this  study. 
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Figure  15.  Model  Sensitivity  to  S*  Variable  Changes  vs  T:  Model  results  when 
F  and  S'  are  held  constant,  but  the  variables  within  S’  are  varied. 


Here,  the  freezing  efficiency,  £*,  is  constant  as  are  the  other  parameters, 


while  A T  has  been  changed. 
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Figure  16.  Model  Sensitivity  to  S*  Variable  Changes  vs  S:  Model  results  when 
F  and  S'  are  held  constant,  but  the  variables  within  S'  are  varied. 
Here,  the  freezing  efficiency,  £*,  is  constant  as  are  the  other  parameters, 
while  6S  has  been  changed. 
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Table  7.  VALUES  FOR  S*  SENSITIVITY  MODEL  RUNS 


Variables 

Value  Pairs  Used 

A  T 

-.287,  -.574,  -1.15,  -1.72,  -2.01,  -2.30,  -2.58,  -2.87,  -3.44,  -4.02,  -4.59,  -5.16, 

-5.74 

m 

-0.05,  -0.10,  -0.20,  -0.30,  -0.35,  -0.40,  -0.45,  -0.50,  -0.60,  -0.70,  -0.80,  -0.90, 

-1.00 

C.  MODEL  SENSITIVITY  TO  NON-DIMENSIONAL  PARAMETERS 

To  determine  if  most  of  the  model  variability  could  be  accounted  for  by  the  three 
non-dimensional  parameters,  model  runs  were  performed  where  the  resultant 
non-dimensional  parameters  were  not  changed,  but  the  individual  variables  in  the 
parameters  were  varied. 

1.  Model  Sensitivity  to  S*  Variable  Changes 

For  this  study  the  two  model  variables  comprising  S',  AT  and  AS,  were  varied 
so  as  to  keep  the  resultant  value  of  S'  constant.  Table  7  lists  the  values  used  for  AT 
and  AS  for  these  model  runs.  The  values  for  all  the  other  variables  and  initial  conditions 
are  as  in  listed  in  Table  5  under  normal  values. 

Figure  15  and  16  show  the  results  of  these  runs  where  the  model  variables  are 
kept  constant  except  for  varying  only  the  iwo  variables  which  make  up  S'.  Figure  15 
shows  the  results  relative  to  the  various  values  of  AT.  There  is  a  near  constant  value  for 
F,  S’,  and  F  as  the  values  of  A  T  varied.  Figure  16  show  the  results  relative  to  the 
various  values  of  A  S.  Again  there  are  near  constant  values  for  F,  S’,  and  F  as  AS  is 
varied.  For  example,  in  a  real  situation  the  expected  maximum  freezing  rate  would  be 
the  same  for  the  conditions  of  T  =  —  1.8687,  S  =  34.5,  T0  =  —1.582  ,  S0  =  34.55,  as  well 
as  for  the  other  extreme  where  the  conditions  are  T  —  —1.8687,  S  =  34.5  ,  Tc  =  3.869, 
S0  =  35.5. 
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Net  Surface  Heat  Flux  (Q„) 


(5) 


Model  Sensitivity  to  F*  Variable  Changes  vs  Q:  Model  results  when 
F  and  S'  are  held  constant,  but  two  of  the  variables  within  F,  Q„  and 
h,  are  varied.  Here,  the  freezing  efficiency,  F,  is  constant  as  are  the 
other  parameters,  while  Qn  has  been  changed. 


0.20  0.25  0.30  0.35  0.10  0.15  0.50 


1  > •  I  I  i  l  I 
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Mixed  Layer  Depth  (h)  ( m ) 

Figure  18.  Model  Sensitivity  to  F*  Variable  Changes  vs  h:  Model  results  when 
F  and  S'  are  held  constant,  but  two  of  the  variables  within  F,  Q„  and 
h,  are  varied.  Here,  the  freezing  efficiency,  £\  is  constant  as  are  the 
other  parameters,  while  h  has  been  changed. 
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Mixed  Layer  Depth  (h)  (m) 


Figure  19.  Model  Sensitivity  to  F*  Variable  Changes  vs  b:  Model  results  when 
F  and  S'  are  held  constant,  but  two  of  the  variables  within  F,  h  and 
ul,  are  varied.  Here,  the  freezing  efficiency,  £\  is  constant  as  are  the 
other  parameters,  while  h  has  been  changed. 
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Figure  20.  Model  Sensitivity  to  F*  Variable  Changes  vs  u*:  Model  results  when 
F  and  S’  are  held  constant,  but  two  of  the  variables  within  F,  h  and 
are  varied.  Here,  the  freezing  efficiency,  F,  is  constant  as  are  the 
other  parameters,  while  u?  has  been  changed. 


.  Model  Sensitivity  to  F*  Variable  Changes  vs  u*:  Model  results  when 
F  and  S'  are  held  constant,  but  two  of  the  variables  within  F,  lit  and 
Qn,  are  varied.  Here,  the  freezing  efficiency,  £\  is  constant  as  are  the 
other  parameters,  while  u?  has  been  changed. 
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Figure  22.  Model  Sensitivity  to  F*  Variable  Changes  vs  Q:  Model  results  when 
F  and  S’  are  held  constant,  but  two  of  the  variables  within  F,  id  and 
Q„,  are  varied.  Here,  the  freezing  efficiency,  F,  is  constant  as  are  the 
other  parameters,  while  Qn  has  been  changed. 
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Table  8.  VALUES  FOR  F*  SENSITIVITY  MODEL  RUNS 


Variables 

Changed 

Variabl 

Multiplier  Value  Pairs  Used 

B 

Q„ 

BUB 

0.01,  0.02,  0.05,  0.10,  0.20,  0.50,  1.00,  2.00, 
5.00.  10.0,  20.0,  50.0,  100. 

D 

60  "Y" 

100.,  50.0,  20.0,  10.0,  5.00,  2.00,  1.00,  0.50, 
0.20,  0.10.  0.05,  0.02,  0.01 

h  ,  ul 

D 

60 

0.01,  0.02.  0.05,  0.10,  0.20,  0.50,  1.00,  2.00, 
5.00.  10.0,  20.0.  50.0,  100. 

ul 

0.01,  0.02,  0.05,  0.10,  0.20,  0.50,  1.00,  2.00, 
5.00,  10.0.  20.0,  50.0,  100. 

ul 

0.01,  0.02,  0.05,  0.10,  0.20,  0.50,  1.00,  2.00, 
5.00,  10.0,  20.0,  50.0,  100. 

O' 

0.01,  0.02,  0.05,  0.10,  0.20,  0.50,  1.00,  2.00, 
5.00,  10.0,  20.0.  50.0.  100. 

2.  Model  Sensitivity  to  F*  Variable  Changes 

For  this  study  the  three  model  variables  comprising  F  ,ul,h  and  Qn ,  were  varied 
so  as  to  keep  the  resultant  values  of  F  constant. 

Table  8  lists  the  multiplier  pairs  that  were  used  on  two  of  the  three  variables 
while  the  third  was  kept  constant.  For  each  of  the  model  runs,  the  normal  values  of  two 
of  the  variables  are  multiplied  by  constants  while  the  normal  value  of  the  third  variable 
is  used.  For  example,  the  value  used  for  Q„  is  .01  times  its  normal  value,  the  value  used 
for  h  is  100  times  its  normal  value,  and  the  value  used  for  u?  is  its  normal  value.  The 
resultant  value  of  hQ„  remains  constant  as  Q„  and  h  are  varied. 

Figure  17  and  18  show  the  results  when  only  Q„  and  h  are  varied,  with  ul  held 
constant.  Figure  17  shows  the  results  as  a  function  of  Q„.  Figure  18  shows  the  results 
as  a  function  of  h.  Both  figures  show  no  change  in  F  as  the  variables  Q„  and  h  are 
varied. 
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Figure  23.  The  Limit  of  Freezing:  The  solution  dividing  freezing  and 
non-freezing  regimes  in  terms  of  S’  and  F  The  freezing  regime  is  the 
hatched  area  close  to  the  origin.  In  this  area  the  mixed  layer  in  capable 
of  freezing  with  time.  Outside  it  is  not  possible  for  the  mixed  layer  to 
freeze  until  some  other  forcing  or  stability  parameters  are  changed. 
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Figure  19  and  20  show  the  results  when  only  h  and  ul  are  varied  and  Q„  is  held 
constant.  Figure  19  shows  the  results  as  a  function  of  the  values  of  h.  Figure  20  shows 
the  results  as  a  function  of  the  values  of  ul  .  Both  figures  show  no  change  in  £*  as  the 
variables  h  and  ul  are  varied. 

Figure  21  and  22  show  the  results  when  only  ul  and  Q„  are  varied  while  h  is  held 
constant.  Again  the  results  show  no  variation  of  B  through  the  range  of  values  used. 

These  results  demonstrate  the  value  of  using  the  three  non-dimensional 
parameters  to  study  the  result  of  model  runs.  Only  variations  in  these  parameters  will 
be  considered  and  not  how  those  variations  were  obtained.  This  gives  families  of 
solutions  vice  thousands  of  individual  time-dependent  sohr.ions,  simplifying  the 
understanding  of  how  the  individual  variables  affect  any  prediction. 


D.  THE  LIMIT  OF  FREEZING 

The  regime  where  freezing  cannot  occur  can  be  investigated  using  the 
non-dimensional  parameters.  Appendix  A  contains  the  derivation  of  a  relationship 
separating  the  freezing  non-freezing  regimes: 

S*  =  - L -  (4.4) 

1  +  CjF  +  C2 

This  is  a  relationship  of  only  four  parameters:  the  forcing  parameter,  the  stability 
parameter,  and  the  two  model  constants.  The  reliability  of  the  model  is  dependent  on 
realistic  model  constants,  C,  and  C2.  These  two  model  constants  represent  the  relative 
importance  of  the  two  factors  used  in  computing  the  entrainment  velocity,  IV„  equation 
(3.14). 

Figure  23  shows  the  delineation  between  the  freezing  domain,  the  hatched  area  in 
the  figure,  and  the  domain  of  no  freezing,  the  white  area  in  the  figure,  as  functions  of 
F  and  S'.  From  the  figure  there  is  a  definite  maximum  value  of  S'  where  freezing  can 
occur.  Any  higher  value  of  S',  a  regime  where  the  Ab  is  approaching  0,  will  result  in 
warming  of  the  mixed  layer.  Equation  (4.4)  also  indicates  that  for  any  positive  value  of 
F  there  is  a  S'  where  freezing  can  occur.  For  high  values  of  F  this  equates  to  situations 
where  there  are  very  large  salinity  jumps.  Since  the  coefficient  of  salinity  contraction  is 
so  much  larger  than  the  coefficient  of  thermal  expansion,  this  regime  is 
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Figure  24.  Freezing  Efficiency  as  a  function  of  Forcing  and  Stability:  The 

graphical  depiction  of  equation  (4.7)  where  the  maximum  freezing, 
related  to  the  Freezing  Efficiency,  E\  is  depicted  as  a  surface  above  the 
FxS'  plane.  For  various  conditions  of  forcing  terms  which  show  up  in 
the  F  and  F  terms,  and  stability  terms,  accounted  for  in  the  5'  term, 
a  maximum  freezing  rate  can  be  determined  from  the  F  term. 


not  rare,  and  it  occurs  frequently  is  areas  where  there  has  been  fresh  water  mixed  into 
the  mixed  layer,  i.e.  the  northern  coasts  of  northern  hemisphere  continents  where  there 
is  river  run  off  and  seasonal  ice  pack  melting. 


E.  FREEZING  EFFICIENCY  RELATIONSHIP  TO  FORCING  AND  STABILITY 

A  relation  can  be  obtained  between  the  non-dimensional  parameter  £\  the  forcing 
parameter,  F  ,  the  stability  parameter,  S',  the  model  constants,  C,  and  C2,  and  the 
estimated  amount  of  salt,  or  brine,  extracted  by  freezing.  Appendix  B  derives  a 
relationship  between  the  freezing  rate,  F,  and  other  model  variables  and  constants  and 
is: 


F 


+  Q 


cp 

(S  -  St)  -f- 

Lf 


(4-5) 


u. 


Appendix  C  derived  an  expression  for  —  in  terms  of  other  model  variables  and  is: 

h 


3 

u. 


On 

pLf 


F  g 


+  /?(S-S,) 

S  \  P 


'-p  A 


-f-7T(S-SMCxF  +  C2)  - 


S  -  1 


(4.6) 


Using  equations  (4.5)  and  (4.6),  Appendix  D  derived  an  expression  for  F 
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•D 


Figure  25. 


Freezing  Efficiency  for  Cl  =  1.  and  C2  =  .4:  The  graphical  solution 
for  Freezing  Efficiency  with  the  model  constants  changed  to  C,  =  1.0 
and  C2  =  0.4  . 
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Freezing  Efficiency  for  Cl  =  .5  and  C2  =  .8:  The  graphical  solution 
for  Freezing  Efficiency  with  the  model  constants  changed  to  C,  =  0.5 
and  C:  =  0.8  . 


Figure  27.  Freezing  Efficiency  for  Cl  =  .2  and  C2  =  2.:  The  graphical  solution 
for  Freezing  Efficiency  with  the  model  constants  changed  to  C,  =  0.2 
and  C2  =  2.0  . 
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.  S*2(F*(C,C2L5  -  C,)  +  (1  +  C2){C2LB  -  1))  +  S*(f’C,  +  2  +C2-  C2LB)  -  1 

S*2((l  -  C2LB)(FC}LB+C2LB-  1))  4-  S*(2  —  C}LBF*  —  2C2LB)  —  1  ^ 

where 


•  B  =  /?(S  -  S,) 

Figure  24  is  a  graphical  representation  of  equation  (4.7)  in  S'  x  F  x  F  space. 
This  figure  represents  the  maximum  freezing  that  can  be  expected  from  any  condition 
defined  in  the  F  x  S’  plane. 

From  Figure  24,  when  S'  is  near  zero,  there  can  be  freezing  for  any  situation  where 
there  is  a  heat  flux  up  out  of  the  surface  of  the  mixed  layer.  Also,  for  lower  values  of 
F  a  wider  range  of  salinity; temperature  combinations  will  be  able  to  produce  sea  ice. 

When  the  S'  is  higher  than  .84,  there  is  no  regime  where  freezing  can  occur  because 
there  is  relatively  no  stability/buoyancy  left  between  the  two  layers  and  complete 
overturing  of  the  water  column  will  occur  first. 

Figures  25,  26,  and  27  show  the  effect  changing  the  model  constants  ,  C,  and  C2,  has 
on  the  mode!  output.  From  equation  (3.14),  as  C,  decreases  in  significance  relative  to 
C2  ,  the  effect  of  TKE  generation  from  wind  stress  decreases  and  the  magnitude  of  F 
becomes  less  significant.  This  can  be  seen  in  equation  (E.7)  where  every  F  is  multiplied 
by  a  C,  term.  As  the  value  of  C,  decreases  the  effect  F  has  on  F  vanishes.  There  is  only 
a  slight  change  in  F  as  F  increases  in  Figure  27.  It  does  appear  that  the  volume  under 
the  surface  has  remained  near  constant,  as  the  decrease  in  F  near  0  has  resulted  in 
an  increase  in  F  at  higher  values  of  F. 
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V.  SUMMARY 


A.  CONCLUSIONS 

This  study  has  investigated  the  relationships  of  the  variables  which  define  and  affect 
mixed  layer  dynamics  and  thermodynamics  with  regard  to  the  onset  of  freezing  of  sea 
ice.  The  mixed  layer  and  underlying  water  column  is  modeled  as  a  simple 
one-dimensional  two  layer  system.  The  salient  time  dependent  variables  are  mixed  layer 
temperature,  mixed  layer  salinity,  mixed  layer  depth,  and  mixed  layer  buoyancy. 
Functions  that  relate  the  forcing  terms  and  boundary  values  are  the  entrainment  velocity 
function,  and  the  freezing  rate  function. 

Model  runs  studied  the  sensitivity  of  the  model  to  varying  only  one  parameter. 
These  parameters  included  the  wind  stress,  the  surface  heat  flux,  the  temperature  or 
salinity  jump  between  the  mixed  layer  and  the  deep  layer,  and  the  mixed  layer  depth. 
These  studies  showed: 

•  Increasing  only  the  upward  surface  heat  flux,  Q„,  resulted  in  a  nearly  linear  increase 
in  freezing  rate. 

•  Increasing  the  wind  stress  term,  uf,  resulted  in  a  nearly  linear  decrease  in  freezing 
rate  until  only  a  net  warming  of  the  mixed  layer  was  possible. 

•  Increasing  the  temperature  jump,  A  7",  resulted  in  a  non-lmear  decrease  in  freezing 
rate  until  only  a  net  warming  of  the  mixed  layer  was  poss  Me. 

•  Increasing  the  salinity  jump,  AS,  resulted  in  an  initial  rapid  increase  in  freezing  rate 
followed  by  a  very  gradual  increase  converging  toward  the  maximum  freezing  rate 
expected  for  a  given  surface  heat  flux. 

•  Increasing  the  mixed  layer  depth,  h,  resulted  in  an  initial  rapid  increase  in  freezing 
rate  followed  by  a  very  gradual  increase  converging  toward  a  point  less  than  the 
maximum  freezing  rate  that  would  be  expected  for  a  given  surface  heat  flux. 

Non-dimensional  ratios  formed  from  the  initial  and  boundary  conditions  were 
derived.  A  non-dimensional  parameter,  F,  is  a  ratio  of  the  wind  stress  forcing  of  the 
mixed  layer  of  depth  h  to  the  surface  buoyancy  flux.,  A  parameter,  S',  is  a  ratio  of  the 
buoyancy  due  to  the  temperature  jump  to  the  buoyancy  due  to  the  salinty  jump,  and  is 
a  measure  of  stability.  A  parameter,  £\  is  a  ratio  of  the  freezing  rate  to  the  maximum 
freezing  rate  possible  for  a  given  surface  heat  flux.  The  results  focus  on  the  dependence 
of  F  upon  F  and  S'. 


Model  solutions  revealed  no  sensitivity  to  changes  in  the  model  variables 
constituting  the  non-dimensional  parameters,  if  the  value  of  the  non-^'mensional 
parameters  is  held  constant.  This  indicates  that  the  model  may  be  completely 
determined  by  the  non-dimensional  parameters  in  certain  domains. 

A  relationship  between  the  forcing  parameter,  F  ,  and  the  stability  non-dimensinal 
parameter,  S‘,  was  derived  which  defined  the  oceanic  domains  where  freezing  could 
occur  and  where  freezing  could  not  occur. 

An  analytic  solution  of  the  model  was  derived  for  the  case  of  freezing  conditions. 
This  solution  showed  the  freezing  efficiency  non-dimensional  parameter,  F,  was  a 
function  of  F,  S',  physical  constants,  and  two  model  constants  for  entrainment  rate. 
This  solution  was  then  displayed  and  analyzed  for  significant  features,  showing: 

•  The  freezing  efficiency  is  highly  dependent  on  the  value  of  S‘  ,  the  stability 
parameter. 

•  If  the  stability  parameter  is  large,  there  can  be  no  freezing. 

•  If  the  stability  parameter  is  near  zero,  there  will  always  be  a  situation  where 
freezing  can  occur  regardless  of  the  value  of  the  forcing  parameter. 


B.  RECOMMENDATIONS 

The  model,  the  results  of  this  study,  and  the  analytic  solutions  can  be  valuable  in 
understanding  the  interactions  between  the  variables  comprising  mixed  layer 
thermodynamics  and  dynamics,  and  their  relationship  to  the  freezing  process.  These 
results  could  be  used  further  to  study  the  following: 

•  The  possible  cvciic  effect  of  cooling  from  above  and  heating  from  below  as  studied 
by  VVelander  (1977). 

•  The  negative  feedback  effect  of  the  formation  of  ice  on  the  freezing  rate  as 
discussed  by  Chu  and  Garwood  (1988). 

•  The  freezing  rate  as  the  drag  coefficient  decreases  from  the  values  of  the  open 
ocean  with  the  initial  formation  of  ice,  and  then  increases  with  the  growth  of  the 
ice,  based  on  the  results  of  Guest  and  Davidson  (1991). 


76 


APPENDIX  A.  THE  LIMIT  OF  FREEZING 


The  transition  region  between  freezing  and  not  freezing  will  be  investigated.  In  this 
regime,  the  mixed  layer  temperature  is  at  the  freezing  temperature  but  there  is  no 
freezing.  In  this  regime,  mixed  layer  temperature  cannot  change  in  time.  This  will 
separate  the  regimes  where  the  mixed  layer  f’ools  with  time  and  where  the  mixed  layer 
warms  with  time.  For  this  investigation  evaporation,  precipitation,  melting,  and  freezing 
are  all  assumed  zero. 


A.  THE  MODEL  EQUATIONS 

For  the  regime  of  no  change  of  mixed  layer  temperature  with  time, 

^  -  0 
Cl 


the  point  where  freezing  rate  is  zero, 

F  =  0 

and  there  is  no  melting,  evaporation,  or  precipitation,  giving 

E=F=P=M=  0  (A.l) 

and  when  the  mixed  layer  is  at  the  freezing  temperature, 

T  =  7) 

the  model  equations  will  be  investigated  for  a  solution. 

For  the  condition  of  no  change  of  temperature  with  time,  the  temperature  flux  at  the 
base  of  the  mixed  layer  and  the  temperature  flux  at  the  top  of  the  mixed  layer  must  be 
equal, 

7%7(0)  «  7V(-/i) 

From  equation  (3.10),  the  thermal  flux  at  the  base  of  the  mixed  layer  is 
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T  w'(  —hi)  =  -  A  T\Ve 


and  from  equation  (3.7)  the  thermal  flux  at  the  top  of  the  mixed  layer  is 


7V(  0)  =  - 


This  gives 


=  -  A  T 


Rearranging  equation  (A. 2)  for  W,  gives, 


W  =  — — — 

'  pcM 


Equation  (3.14)  with  equation  (A.l)  gives 


/  3 

-  is-  c~ 


-  C2ag- 


Equating  equations  (A. 3)  and  (A.4)  together  gives 


Qn  _  1  ( r  r  Qn 

pc„&T  -  Ab  iCi  h  “  PCp 


Equation  (A. 5)  can  be  manipulated  to  give 


A  b  = 


AT  pc. 


:-(c4 


r  Ml 
C2aS  pc. 


—  AT  I  C, 


-  C2<xg 
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Ab  =  agAT 


From  equations  (3.9)  and  (A.l) 


" 

_ 

“ 

3 

c, 

U. 

-  c2 

hagQ" 

L  pcp  J 

i'^(0)  =  - 


Substituting  equation  (A. 7)  into  equation  (A.6)  gives 


Ab  —  <xgA  T 


3 

u. 


-hb'w'{0) 


-  C, 


Using  the  definition  of  F,  equation  (4.2),  equation  (A. 8)  becomes 

Ab  =  agAT(-  C,F*  -  C2) 

Using  equation  (3.13)  equation  (A. 9)  becomes 

agAT  -  pgAS  =  agAT ( -  C,  F  -  C2) 

Rearranging  equation  (A.  10) 

agAT(l  +  C,F*  +  C2)  =  /?gAS 


(A.6) 


(A. 7) 


(A.8) 


(A. 9) 


(A.  10) 


or 


1  +  C,F  +  C-)  = 


agAt 


(A.l  1) 


B.  THE  FREEZING  LIMIT  SOLUTION 

From  the  definition  of  S’,  equation  (4.3),  equation  (A.l  1 )  can  be  rewritten 

1  +  C,F  +  C2  =  -T-  (A. 12) 

S 

Solving  equation  (A.  12)  for  S' 
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(A.  13) 


1  +  CjF  +  C2 

This  gives  a  relationship  of  four  parameters:  the  forcing  parameter,  the  stability 
parameter,  and  the  two  model  constants.  This  indicates  that  the  reliability  of  the  model 
is  dependent  on  realistic  model  constants. 

Equation  (A.  12)  can  be  solved  for  F 

f  =  -  (1  +  C2))  (A- 14) 
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APPENDIX  B.  AN  EXPRESSION  FOR  THE  FREEZING  RATE 


A.  THE  CONDITIONS 

It  may  be  possible  to  derive  an  analytic  solution  for  the  freezing  rate  from  the 
equations  used  with  the  model.  This  solution  will  be  useful  for  comparing  model  results 
with  analytic  results,  but  also  useful  for  clarifying  relationships  between  the  model 
variables. 

The  main  purpose  for  finding  an  analytic  solution  would  be  to  find  regimes  where 
we  would  expect  freezing  or  regimes  where  we  would  not  expect  freezing.  This  has  very- 
useful  applications  in  ocean  areas  where  preconditioning  may  or  may  not  permit  rapid 
freezing  to  occur.  Naval  operations  may  be  extended  if  conditions  are  not  met  for 
freezing  when  freezing  is  a  limitation  to  operations,  or  operations  may  need  to  be 
curtailed  early  due  to  forecast  conditions  which  could  lead  to  freezing  conditions. 


B.  ESSENTIAL  EQUATIONS 

The  basic  equations  which  will  be  needed  are  the  equation  defining  freezing  rate, 
entrainment  velocity,  and  the  equations  defining  the  non-dimensional  parameters. 
Equations  (3.18)  and  (3.14)  will  be  used. 


C.  THE  DERIVATION 


From  equation  (3.18),  F,  the  Freezing  Rate,  is  defined  as 


F 


+  AT 


We 


From  equation  (3.14),  IV„  the  Entrainment  Velocity,  is  defined  as 
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XV. 


Ab 


3 

U* 


h 


PgF(S  -  St) 


Substituting  equation  (3.14)  into  (3.18)  for  IV,  gives 


F  = 


Qn_ 

pcp 


+  AT\ 


1 

Ab 


PgF{S  —  Sj) 


or 


r-  Qn  r  AT  “•  CP 

F  =  -—  +  C, 


pLf 

Collecting  terms  gives 


Ab  h  L 


'/ 


^  AT  Qn  ^  AT  cp 

c2ag~nr~rr  +  c2figF[s-s,)- 


Ab  pLf 


Ab  Lf 


From  equation  (3.13)  and  (4.3) 


Ab  =  erg  AT  -  P  g  AS 
AT 


AT 


or 


Ab  r,  i  l 

Tf  -  ■*[>-  y\ 

Multiplying  both  sides  of  equation  (B.2)  by  (B.3)  gives 


(B.3) 


1 

* 

1 _ 

«g 

L  J 

-  C2pg{S  -  S.) 


Qn 

S  -  I 

+  C2crg 

pLj 

ag 

S. 

cn  u, 

+  c'i~ 


(B.4) 
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or 


F\ 


C2 


£ 

a 


+  C, 


cp  ul 

LjOLg  h 


(B.5) 


Finally 


F 


+  C2 


+  c, 


(B.6) 


D.  POINTS  OF  INTEREST 


From  equation  (B.6),  as  S'=s>0  then 


or  from  equation  (4.1) 


£*  =>  1 


For  F  =  0 


Qn_  (S'  -  1 
pLf  y  S. 


C, 


or 


5  -  1 
5. 


+  C2  -  C, 


3 

u* 


(B.7) 


P  If 
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From  equation  (3.9),  when  F  =  0,  then 


b'w’{  0) 


ag 


Qn_ 

PCp 


Substituting  (B.8)  into  equation  (4.2)  gives 


* 

F 


3 

U« 


*g 


Oil 

pcp 


Substituting  equation  (B.9)  into  (B.7)  gives 


C,(-  F) 


Solving  for  S'  gives 

S’  -  ( - i-^) 

v  1  +  C2  +  C,  F  / 

This  is  the  same  result  as  equation  (A.  13),  the  distinction  between  regimes 
and  not  freezing. 


(B.8) 

(B.9) 

(BIO) 

(B. 11) 
of  freezing 
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APPENDIX  C.  AN  EXPRESSION  FOR  U* 


A.  BEGINNING  EQUATIONS 


An  expression  for  u }  can  be  obtained  from  the  model  equations  and  from  the 
expression  of  F.  From  equation  (4.2) 


hb'w'i  0) 


From  equation  (3.9) 


b'w'{  0) 


-  *g 


Qn_ 

pcp 


+  0gF(S  -  S,) 


From  equation  (B.6) 


F 


Qn 

pLf 


S  ~  1 
S. 


+  c. 


„  CP  «* 

+  C,  — -  ~ T~ 

LfCtg  h 


B.  THE  DERIVATION 


Equation  (4.2)  can  be  rewritten 


Fm  b'w'(0) 


3 

U. 

~h 


and  substituting  equations  (3.9)  and  (B.6)  in  gives 
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«S-y-  +  MS -Si) 


Qn 

pLf 


S  -  1 


+  c2  ]  +  c, 


3 

u. 


txgLf  h 


S  -  1 


A  '-p 


-  -  C2T-f(S-S,) 

s  h 


3 

U» 


Multiplying  both  sides  by  the  denominator  gives 


Qn  (  S'-  1  r 

pc7  -7  C2  «  U{S  Sl) 


+  F 


MS  - 


s'  —  1  c\  ,  c 
_  +  c2j  +  ct  agLf  h 


Qn 


(C.l) 


w.  /  s*  -  1  c  P  s) 
=  —  -c=“  l/s  s‘\ 


or 


»CP  s' 


pLf 


+  F 


Ms-s')-$’('STL  +  Cl)  +  C'^^T{S-Sl) 


S  -  1 

h  s' 


6  cp  u, 

cA-r:^-^ir 


Gathering  terms 
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C.  POINTS  OF  INTEREST 

From  equation  (C.2)  some  points  can  be  drawn.  First,  when  S'=>0  then 

Qn  •  1" 

-  -  -  + 

This  means  that  when  the  mixed  layer  has  no  temperature  jump  between  it  and  the 
mixed  layer,  i.e.  AT  =  0,  then  the  freezing  rate  is  only  a  function  of  the  thermal  heat  flux 
out  of  the  mixed  layer.  The  reduction  of  TKE  is  accomplished  only  thru  the  entrainment 
of  more  dense  saline  deep  water. 
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u . 

When  S'=s-1,  —  =>0  or  h=>o o  This  means  that  any  buoyancy  flux  at  the  surface 
cannot  be  countered  by  a  buoyancy  flux  at  the  bottom  of  the  mixed  layer  and  the 
subsequent  condition  would  be  an  over  turning  of  the  water  column. 
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APPENDIX  D.  FREEZING  EFFICIENCY 

A.  THE  EQUATIONS 

A  expression  for  E  will  be  derived  from  equations  (B.6)  and  (C.2).  The  resulting 
expression  will  be  a  function  of  the  model  constants,  the  salinity,  and  the 
non-dimensional  parameters  F  and  S’.  This  will  allow  an  analysis  of  the  possible 
maximum  freezing  rate  possible  for  any  initial  conditions  and  any  boundary-  conditions. 


B.  DERIVATION 

Equation  (C.6)  gives  the  expression  for  F  as 


__  Qn_ 

r  S'.  -  1  ,  r  1 

„  cp  M* 

pLf 

s.  + 

1  Ffixg  h 

\  1 

1  * 

V) 

l _ 1 

B  Cp 

-  Cj  (S  -  S,)  -f- 

or 


F  = 


Qn_ 

plf 


(s'- l  +  c3s-)  +  c, 

.  B  .  Cp 

S  -  1  -  C2  £  (S  -  Si)  s  -Z- 


(D.l) 


Equation  (C.2)  gives  the  expression  for  -7-  as 

n 


3 

u. 


Qn_ 

pLf 


Fg 


S  -  1 


+  P(S~S,) 


v  P 


—  -(S-S^QF  +  C2)  - 


s  - 
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or 


(s'->)(«77  +  «s-sj) 

~-y(S-S,)S'(C,F  +  C-j  -  s'  +  I 
Lf  ' 

To  simplify  the  writing  of  the  equations,  let 


and 

B  =  p(S-Si) 

Equation  (D.  1 )  becomes 

/  *  \  7  *  Urn 

-  +  +  c<  g  s  IT 

f  —  - 1 - 

S*  -  1  -  c2lbs 

Equation  (D.2)  becomes 

LBS‘{C\F  +  C2)  -  S*  -  1 
Substituting  equation  (D.4)  into  (D.3)  gives 

r  _  _ 

5*(1  -  C2LB)  -  1 


(5*  —  IVy-  +  B) 
S’iLBC^F  +  C2-  1)+  1 


Um  -  r* 
h  pLf  g 


(D.2) 


(D.3) 


(D.4) 
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or 


(5  (1  +  C2)-  1)  -  C,  S  F 


Qn 

pLf 


(S  -  1)(1  +  LB) 


S  (L5(C,£  +  C2)  -  1)  +  1 


S  (l  ~  C2LB)  -  I 


<D.5) 


From  equation  (4.1) 


E  = 


-Qn 

pLf 


Equation  (D.5)  can  be  written 


(S*(l  +  C2)-l)  -  C,5*E’ 


£  = 


(S  -  1)(1  +  LB) 


S  (LB(CXF  +  C2)  -  1)  +  1 


5(1  —  C2  L  B)  —  ] 


or 


(S’(l  +  C2)  -  l)(S*(£fi(Cj£  +  C2)  -  1)  +  l)  -  CjS*£*(5*  -  1)(1  +  L  B) 
(S*(l  -  C2  LB)  -  l)(5’(L£(C,£*  4-  C2)  -  1)  +  l) 


or 


£  = 


S*2(l  +  C2)(L£(C,£  +  C2)  -  1)  +  5*((1  4-  C2)  —  ( LB(C}F  +  C2)  -  l))  -  1 

S*2((l  -  C2LB)(LB(CXF ’  4-  C2)  -  1))  +  S*((l  -  C2LB)  -  (££(C,£*  +  C2)  -  1))  -  1 

(D.6) 

CjS*£*(S*(l  +  LB)  -  (1  4-  LB)) 

5*2((1  -  C2LB){LB{CXF *  +  C2)  -  1))  +  S'((l  -  C2LB)  -  (££(C)£*  1-  C2)  -  1))  -  1 
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Finally,  equation  (D.6)  can  be  written 

.  S,2(F*(C,  C2LB  -  C|)  +  (1  +  C2)(C2LB  -  1))  +  S'(FC}  +2  +  C2-  C2LB)  -  1 

S*2((l  -  C2LB){F’CxLB  +  C2LB  -  1))  +  S\ 2  -  CXLBF  -  2 C2LB)  -  1 


<D.7) 


where 


•  F  = 


•  F  = 


•  S'  = 


pLF 

-2. 

M? 


hb'w'{  0) 

agAT 

AgAS 


•  L  - 


aLf 


B  =  P(S  -  St) 


•  C|  and  C2  are  model  constants 
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